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A&s*rac<-Because material junctions are commonplace on structures whose radar cross section is 
o interest it is essential that their scattering properties be adequately characterized. The stan- 
dard impedance boundary condition (SIBC) has been employed in the past along with function 
theoretic techniques to develop simple scattering models of material junctions with thin and/or 
high loss slabs. To extend these models to more general slabs, generalized impedance boundary 
conditions (GIBCs) and generalized sheet transition conditions (GSTCs) have been proposed 
Unfortunately, the solutions obtained with these are usually non-unique in the form of unknown 
constants and although the constants have been resolved for a few special cases, previous efforts 
were unable to determine them in the general case. 

This report examines the problem of the plane wave diffraction by an arbitrary symmetric 
tW(>dimensional junction , where Generalized Impedance Boundary Conditions (GIBCs) and Gen- 
erahzed Sheet TVans.tion Conditions (GSTCs) are employed to simulate the slabs. In chapter 2, 

re uffL a GIBC JcSTC £" mU ' tllayer P lanar slabs * arbitrary thickness and the 

S C Rrf/rSr T C u ' T ‘° e ® C,ents are com P ared with e *act counterparts to evaluate 
the GIBCs/GSTCs. In chapter 3 the plane wave diffraction by a multilayer material slab recessed 

MatrivP* y i PlanC 18 formulated and solv « d via the Generalized Scattering 

. ormu a ion (GSMF) in conjunction with the dual integral equation approach. Various 
scattering patterns are computed and validated with exact results where possible. 

rni V A T 4 ’ ‘ h ® diffraction by a mate bal discontinuity in a thick dielectric/ferrite slab is 
considered by modelling the constituent slabs with GSTCs. A non-unique solution in terms of 
un nown constants is obtained, and these constants are evaluated for the recessed slab geometry 
of chapter 3 by comparison with the solution obtained therein. Several other simplified case's 
f a :;r prese " e an discussed. In chapter 5 an eigenfunction expansion method is introduced 
o determine the unknown solution constants in the general case. This procedure is applied to 

' T C f. apter4 ’ and fettering patterns are presented for various slab junctions and 

compared with alternative results where possible. Chapter six presents a short summary of this 

report and some recommendations for future work. 



TABLE OF CONTENTS 


LIST OF FIGURES 

LIST OF TABLES 

CHAPTER 

I. INTRODUCTION . 


iii 

viii 

1 .. 


II. DERIVATION OF GENERALIZED TRANSITION/ 
BOUNDARY CONDITIONS FOR PLANAR MULTIPLE 
LAYER STRUCTURES 

2.1 Derivation of Infinite Order Conditions 

2.2 Low and High Contrast Approximations for Matrix Element 

Operators 

2.3 Evaluation of the Boundary/Transition Conditions 

2.4 Concluding Remarks 

III. DIFFRACTION BY A MULTILAYER SLAB RECESSED 
IN A GROUND PLANE VIA GENERALIZED IMPEDANCE 
BOUNDARY CONDITIONS 

3.1 Description of GSMF Procedure 

3.2 Plane Wave Diffraction and Mode Coupling 

3.3 Reflection and Launching of a Waveguide Mode 

3.4 Computation of Spectra for Material Insert in a Perfectly Con- 
ducting Ground Plane 

3.5 Specialization to the GIBC Representation 

3.6 Numerical Results 

3.7 Summary 


14 

22 

27 

43 


i 


n . general solution of the diffraction by a 
materia I discontinuity in a thick dielectric, 

FERRITE SLAB 


4.1 Dual Integral Equation Formulation 7 (; 

4.J Solution ot the Dual Integral Equations ,v> 

4.4 Discussion of the Solution and Some Applications ,s<; 

4.3.1 Diffraction by thin single layer discontinuous slabs . 8b 

4.3.2 Diffraction by a resistive- resistive junction ;)•> 

4.3.3 Diffraction by grounded metal-dielectric junctions . 03 

4.3.4 Diffraction by a thin dielectic/ferrite half plane . . . DO 

4.4 Conclusion , in 


V. RESOLUTION OF NON-UNIQUENESS ASSOCIATED WITH 
THE GIBC/GSTC SOLUTION H 


5.1 Modal Decomposition for the Symmetric Slab Problem .... 103 

5.1.1 Exterior Cross Section Functions l 05 

5.1.2 Interior Cross Section Functions }()7 

5.2 Recasting of the Dual Integral Equation Solution for a Mate- 
rial Discontinuity i ii 

5.3 Determination of the Constants p>-> 

5.4 Validation of the Solution i.>o 

5.5 Summary 


VI. SUMMARY 

APPENDICES . . . 
BIBLIOGRAPHY . 


ii 



LIST OF FIGURES 


Figure 

1.1 Structures modelled by SIBCs - 

1.2 Resistive/conductive sheet simulation of thin slab 4 

2.1 Infinite multilayer slab 1 . 1 

2.2 (a) Multilayer coating on ground plane, (b) Equivalent opaque sheet. 12 

2.3 (a) Zero thickness resistive and conductive sheet similation of multi- 
layer slab, (b) Distributed resistive and conductive sheet simulation 

of multilayer slab 13 

2.4 Coordinate system for reflection and transmission coefficient deriva- 
tion 27 

2.5 E y polarization low contrast simulation of a single layer having rel- 
ative constitutive parameters t = 2, f.t = 1.2 30 

2.6 E y polarization low contrast simulation of a single layer having rel- 
ative constitutive parameters t = 3.5, p* = 2 31 

2.7 E y polarization high contrast simulation of a single layer having 

t = 3.5, /i = 2 32 

2.3 E y polarization high contrast simulation of a single layer having 

e = 5, ji = 3 33 

2.9 H y polarization low contrast simulation of a single layer having rel- 
ative constitutive parameters e = 2, /j, = 1.2 36 

2.10 H y polarization low contrast simulation of a single layer having rel- 
ative constitutive parameters t = 3.5, /z = 2 37 

2.11 H y polarization high contrast simulation of a single layer having 

e - 3.5. fi = 2 33 

iii 


2.12 

2.1:5 

2.14 


2.15 


H y polarization high 
t: = 5. = 3 


contrast simulation of a single layer having 


Simulation of a three layer slab using a 14“ order condition. 

■Simulation of a three layer coating on a PEC ground plane usin 
14 order condition 


a 


Simulation of a three layer coating on a PMC ground 
14 l order condition 


plane using 


cl 


40 


41 


3.1 

3.2 

3.3 


3.4 


3.5 


3.6 


3.7 


(a) Multilayer slab recessed in a PEC ground plane, (b) Represen- 
tation of slab as surface with reflection coefficient R ‘ 


Illustration ot recessed stub geometry (aj and associated subprob- 
uns. (b) direct diffraction, (c) mode coupling, (d) stub reflection, 
( ) m ° de reilectlon at the waveguide mouth, (f) mode launching. 

^polarization backscatter echo width for a material insert with 

f. rv, c “ 2 f • ? ' 0001, ^ = L2 modeled by low contrast GIBCs 
(see fable 3.1 for an explanation of the legend entries), (a) direct 
diffraction component (b) modal component (c) composite. . . 

//^polanzation^backscatter echo width for a material insert with 
,• T 7, 6 .. 7° J - 0001 ’ ^ - 2 modeled by low contrast GIBCs 

,. ff a e - ’ for an explanation of the legend entries), (a) direct 
diffraction component (b) modal component (c) composite 

^polarization^ backscatter echo width for a material insert with 

(see Table 7*? ~ M = 2 modeIed b - v hi S h contrast GIBCs 

( e fable 3._ tor an explanation of the legend entries), (a) direct 

iffi action component (b) modal component (c) composite 

^polarization backscatter echo width for a material insert with 

Lee TaKlVo f ' ~ J ‘ 00 ° 1 1 ’ ** = 7 mod eled by high contrast GIBCs 
(see fable 3.2 for an explanation of the legend entries) 

Hz P° lari ^tion backscatter echo width for a three-laver material 

•H - , W 000l (ri = tV 1 = \\ - jm0h = 7 >’ = -4A, c 2 = 

To 3-0001, - *), (r 3 = 2A, e 3 = 2. - >.0001, fi 3 = 1.2). I n the 

?° e ?. ^ Utl ^ a ’k,c(d,e), a anc * b denote the approximation in of 

the high uoustuust layers 1 and 2, c denotes the approximation in 

r of the low constast layer 3, while d is the order of the composite 

f-im. and e is the total number of modes. 


45 


46 


63 


64 


66 


70 


iv 



ot 


4.1 


(a) Distributee! sheet, (b) Distributed slieet discontinuity, 


4.2 Geometries modeled by a discontinuous distributed sheet, (a) Ma- 
terial half-plane, (b) Material-material join, (c) Grounded join. ... 75 

4.3 (a) Illustration of C contour in the complex a plane and (b) complex 

A plane 79 

4.4 Thin slab structures and their distributed sheet representation (a) 

Material-material junction, (b) Material-metallic join, (c) Material 
half-plane, (d) Equivalent distributed sheet *7 


4.5 (a) Recessed slab (PEC stub) (b) Grounded slab with truncated 

upper plate (c) Recessed slab (PMC stub) (d) GIBC sheet 94 

4.G II. polarization backscatter echo width for a material insert having 
w = .04A. 6 = 2— /'.OOOl, u = 1.2 modeled with O (w) low contrast 
GIBCs 97 - 


4.7 H. polarization backscatter echo width for a material insert having 
w — . 10A, 6 = 2 — >0001. n - 1.2 modeled with 0(w 3 ) low contrast 
GIBCs; bio = — .4578+ j. 2593, b n = — .0408 — >.011 1, b 20 = —.0401 — 
>0388 98 

5.1 Symmetric multilayer slab discontinuity illuminated by a plane wave 102 


5.2 Symmetric multilayer slab with irregular termination illuminated by 

a polarized field excitation 103 

5.3 (a) C contour in the complex A plane, (b) Deformation for region 1 

integrals, (c) Deformation for region 2 integrals 115 

5.4 Indexing scheme for constants 124 


H z polarization backscatter echo width for a recessed slab with r = 

.4, e = 2, n = 1.2 modeled by O(t) second order low contrast GIBC 
(see Table 5.2 for constants) 133 

5.6 II. polarization backscatter echo width for a recessed slab with r = 

.1.6 = 2. /i = 1.2 modeled by 0(t 3 ) fourth order low contrast GIBC 
(see Table 5.2 for constants) 131 


v 


/ P° larization backscatter echo width for a recessed slab with - = 

i<ee Tah 1 3 modelecl b - v 0 ^ _1 ) second order high contrast GIBC 

l^ce table o.J for constants) 


^polanzanon backscatter echo width for a material half-plane with 

't i I ' l v - 1-2 modeled by 0{t) low contrast GSTC (see 
Table o.2 for constants). ' 

//^polarization backscatter echo width for a material half-plane with 

-p , i -Vf ^ — mo< leled by O(r) low contrast GSTC (see 
fable o.2 for constants) ' v ^ 

poUn^tion backscatter echo width for a material half-plane with 

T , ] -'-A ^ 3 ’“ mo ^ ec l by O(r) low contrast GSTC (see 

Table o.2 for constants) 1 


^^polarization backscatter echo width for a material half-plane with 

CSTCI^T l/- = , T 2 m ° deled by ° (r) and °< T ’> low ^trast 
OblL (see Table o .2 for constants) 

//^polarization backscatter echo width for a material half-plane with 

rV-G / °* /Z = 3 modeled b y 0(k°) (even symmetry) and 

•7 °f d } s >7 metr y) high contrast GSTC. The free spa^e side 

is moaeled bv 0 (t) (even svmmptrvl /■)/_ . 2 \ / i i 
_ e a c -rVi symmetry) G>(r odd symmetry) ow 

contrast GSTC. (see Table 5.3 for constants) HO 

//.^polarization backscatter echo width for a material half-plane with 

0 ,~ i, ; A r ° ' " = ) u modeled ^ 0(«°) (even symmetry) and 

is J ft Tf 1 !' 1 h ' gh COntrast GSTC ' Th e free space side 

LTatrCTC r (r ’if' rn” Symmetr - V t 0(r*) (odd symmetry) low 
contrast GSTC. (see Table 5.3 for constants) 141 

H . z P olanz ation backscatter echo width for a two laver .07A svmmet- 
nc materia half-plane with r, = .005, e, = 5, p, = 3 and r, = 03 

, 2 ‘ La 5' er 1 ls modeled by 0(< c°) GSTC, layer 2 and the 
free space layer are modeled by 0(r) (even symmetry') 0(r=) (odd 
symmetry) low contrast GSTC (see Table 5.3 for constants). . . , 42 

//.. polarization backscatter echo width curves for different thickness 
mater, al half-planes with e = 2, p = 1.2. All cases are modeled by 
0( . ) low contrast GSTC (see Table 5.3 for constants) ' m 


vi 



5.16 H z polarization backscatter echo width curves for different thickness 
grounded slabs with e = 5. /.t = 3. The number to the left of the 
colon denotes the O(r) of the l.c. GIBC for the free space side, and 
the number to the right of the colon denotes the 0(k 1 ) ot the h.c. 

slab GIBC 144 

5.17 H. polarization backscatter echo width curves for different thickness 
grounded slabs with t = 11, n = 7. The number to the left of the 
colon denotes the 0{t) of the l.c. GIBC for the free space side, and 
the number to the right of the colon denotes the 0{k *) ot the h.c. 

GIBC for the right hand slab 1 

5. IS H z polarization backscatter echo width curves for different thickness 
grounded junctions, with = 2. /ii = 1.2 and e 2 - 5. f-tj = 3. The 
number to the left of the colon denotes the 0(t) of the l.c. GIBC tor 
the left hand slab and the number to the right of the colon denotes 
the 0 (k - 1 ) of the h.c. GIBC for the right hand slab 1 16 

5.19 H z polarization backscatter echo width curves for different thickness 
grounded junctions, with Ci = 5, Hi — 3 and ( 2 = 11, // 2 = ”• 

The number to the left of the colon denotes the 0 (k _ 1 ) of the h.c. 

GIBC for the left hand slab and the number to the right of the colon 
denotes the 0 (k - 1 ) of the h.c. GIBC for the right hand slab 147 

5.20 H z polarization backscatter echo width curves a grounded junctions, 
with ti = 2, /.ix = 1.2, 7-j = .025 and e 2 = 5, jx 2 = 3. The first number 
to the left of the colon denotes the O(t) of the l.c. GIBC for the 
left hand slab and the first number to the right of the colon denotes 
the 0(/v _1 ) of the h.c. GIBC for the right hand slab. The second 
number on either side of the colon denotes the O(t) of the l.c. GIBC 

for the free space slab needed to give the two sides the same thickness. 14S 

5.21 H , polarization backscatter echo width curves a grounded junctions, 
with ti = 5, m = 3, t 1 = .05 and e 2 = 11, /i 2 = 7- The first number 
to the left of the colon denotes the 0(/c -1 ) of the h.c. GIBC for the 
left hand slab and the first number to the right of the colon denotes 
the 0{k~ 1 ) of the h.c. GIBC for the right hand slab. The second 
number on either side of the colon denotes the 0(t ) of the l.c. GIBC 

for the free space slab needed to give the two sides the same thickness. 149 


list of tables 


e 




^piuAimauons to q tJ Q per 




High contrast approximations to q tJ operators. 

i° itr de r- for * ** < - 

condition, the three numbers ofVh^ j contrast boundary 

order of the approxima.LTn th ^ r T," lnd '“" 
boundary condi, ' M - “ d ^ 

- » rj r g ?f r des for * ^ 
b ftz zz°:z 

^S:r^ttH hw sirai,M '°" ° f ** 

functions Appendtx 1 for definition of M t S p, it 

S r s,raiia, '°" ** 

functions ‘ PPendix 1 for definition of Af + sp ij t 

Asymptotic behavior of integrand fnrtt f 

mte 0 rand for the functions F dtf} {x = 0 ±, y). { 

Values of so, u, .on constants for curves presented in Figures o to 12 . 

Values of so, ion constants for curves presented in Figures 12 to 13 



CHAPTER I 


INTRODUCTION 


The use of non-metallic materials is now commonplace on airborne vehicles, and 
frequently two material slabs of different composition will abut each other to form 
a junction. For example, thin radar-absorber material (RAM) strips of different , 
composition are often joined end-to-end on a metal surface for radar cross section 
reduction. Also, ill microstrip antenna configurations a discontinuity in substrate 
composition is used to suppress unwanted surface waves and in many cases a planar 
array is terminated at a metal-dielectric junction. The electromagnetic effectiveness 
of structures and devices such as RAM coatings or microstrip antennas is influenced 
by the scattering behavior of any material junctions present. It is therefore important 

to obtain a characterization of their behavior. 

The scattering behavior of a material slab junction is revealed by an examination 
of its plane wave diffraction, and this is the overall topic of the dissertation. For 
discussion purposes, it is useful to divide these into non-penetrable and penetrable 
slab junctions. The scattering problem associated with non-penetrable junctions 
was aided greatly by the introduction of standard impedance boundary conditions 
fSIBCs) [35], introduced to model thin metal backed coatings and layers of high 
loss (see Figure 1.1). Conceptually. SIBCs are first order boundary conditions ' 
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Figure 1.1: Structures modelled by SlBC's. 
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discontinuity must remain finite. This two step approach in computing the diffraction 
from slab junctions has been successfully exploited by many researchers [34] [IS] rli ; 
[27; r 30j. However, the approach is predicated on the validity of the SIBC. which is 
restricted to modeling very thin and/or lossy coatings. More general slab junctions 

must therefore be characterized by other methods. 

With regard to junctions formed by penetrable slabs, a close analog to the SIBC 
sheet model is the •resistive'’ and “conductive” sheet simulation [12] [37]. These 
sheets are characterized by simple first order transition conditions \,mch i elate the 
tangential fields across the sheets. In particular, a resistive sheet (see Figuie T-). 
supports an equivalent electric current which produces a discontinuity in the tangen- 
tial magnetic field across the sheet. The proportionality factor relating the tangential 
magnetic field discontinuity to the equivalent electric current is denoted as the ••re- 
sistivity”. Similarly, the conductive sheet model is the dual of the resistive sheet and 
supports an equivalent magnetic current, with the resulting proportionality factor 
denoted as the “conductivity” of the sheet. Like the SIBC models above, resis- 
tive and conductive sheets models of slab half-planes and junctions are amenable 
to Weiner- Hopf methods for the computation of diffraction solutions [1] [44], with 
the same comments given above applying here also. We remark, however, that the 
resistive and conductive sheet models are very restrictive (much more so than SIBC 
models) and cannot be used unless the modeled slabs are very thin and of high pei- 
mittivity and/or permeability. Hence, as in the case of impenetrable slab junctions, 
alternative methods are needed to model more general junctions. 

One such exact approach was employed by Aoki and Uchida [3] to tackle the 
problem of plane wave diffraction from a penetrable single-layer slab junction. Their 
method involved rewriting the junction field components in terms of a Fouiiet se- 
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Figure 1.2: Resistive/conductive sheet simulation of thin slab. 
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nes representation. This approach subsequently led to the generation of Weiner- 
Hopf equations expressed in terms of unknown spectral functions. However, explicit 
expressions for these functions could not be obtained, and the resulting solution 

involved a- cumbersome iterative procedure requiring knowledge of rather complex 
integrals and functions. 

Another possible approach to modeling thicker and more penetrable layered dis- 
continuities is to replace the SIBC and resistive sheet transition conditions (SIC) 
by generalized impedance boundary conditions (GIBC) [15] [17] [23] and generalized 
sheet transition conditions (GSTC) [59] [39], respectively. The GIBC and GSTC are 
equivalent one- and two-sided sheet representations which, unlike their SIBC and 
STC counterparts, display second and possibly higher order derivatives of the field , 
components on the equivalent sheet. These higher order derivatives are responsi- 
ble for the increased accuracy of the GIBC/GSTC relative to the SIBC/STC. and 
the effectiveness of the GIBC/GSTC is roughly proportional to their order. This, 
of course, implies that an adequate GIBC/GSTC modeling of any layered material 
is possible, provided that GIBC or GSTC of sufficient order are employed. Many 
GIBCs and GSTCs have been derived to model all types of single and multiple layers 
[59] [39] [4] [31] [43], However, most of these GIBC/GSTC are either limited in order 
(usually second) or else are only valid for specific geometries. This, of course, points 
to the need for developing more general GIBC and GSTC capable of simulating a 
wide variety of layered slabs and coatings. 

Once a slab junction is represented as a discontinuous GIBC/GSTC sheet, then 
function-theoretic techniques may again be applied to compute the plane wave diffrac- 
tion, as in the SIBC/STC case. However, function theoretic solutions based on 
the application of GIBC/GSTC simulations yield solutions which are non-unique 
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even after the application of the standard edge condition [32] [41] [43] [56], Also 
tec; pi out} is not necessarily satisfied [14] [4] [7], unless this condition is explic- 
itly enforced. Uniqueness is required of any physical solution, whereas reciprocity 
with respect to the transmitter and receiver is necessary when the scattering body 
is electrically passive. As noted in [41]. one may take advantage of this inherent 
non-uniqueness to force a solution which is at least reciprocal, if not unique. The 
non-uniqueness of the solution is usually manifested in terms of unknown constants 
[41] for finite-order GIBC/GSTC and unknown entire functions [32] for infinite-order 
G1BC/GSTC. Most GIBC/GSTC diffraction solutions to date have been obtained 


using second order GIBC/GSTC, and some of these have put forth arguments deal- 
ing with the cancellation of non-physical poles to propose a unique solution [53] [54] - 
[43] [5]. However, these previous efforts offer no method of determining the unknown 
constants for more general GIBC/GSTC simulations. 


If left unresolved, this issue would seriously impede the practical utilization of 
G IBC/GS TC for a characterization of material junctions. In an attempt to resolve it, 
the unknown constant appearing in a second order GIBC solution was recently related 
to the field at the sheet discontinuity [42]. Unfortunately, the edge field is seldom 
known apnori and this relation is not therefore of practical use. Nevertheless, it 


demonstrated that a unique solution may be possible with a GIBC/GSTC simulation. 
An example where it was possible to obtain a unique solution is given by Leppington 
[21], who considered the surface wave reflection by an abrupt change in slab thickness. 
The slab was modelled using second order transition conditions equivalent to those 
given in [59] and [39]. Leppington was able to determine the reflection coefficient 


uniquely in the limiting case of vanishing thickness by matching the interior field far 
from the junction with a static representation of the interior field in the vicinity of 



i 


the junction. This suggests the possibiltiy of working with internal Helds to resolve 
the uniqueness issue, an approach which to date has received little attention and is 
exploited herein. 

The goal of this dissertation is to develop a plane wave diffraction model tor 
general symmetric thick multilayer slab junctions. Four main chapters follow dealing 
with the derivation of the GIBC/GSTC, the formulation and formal solution ot the 
plane wave diffraction by certain GIBC/GSTC approximated slab junctions, and the 
subsequent resolution and explanation of the non-uniqueness phenomenon described 
above. 

In chapter two. arbitrary order GIBC and GSTC are constructed for multilay- 
ered planar slabs of arbitrary thickness. Initially, recurrence relations are derived , 
for the Helds in adjacent layers and are then employed to develop infinite order 
boundarv/transition conditions that are conveniently expressed as a matrix product. 
Approximations to the matrix element operators for low and high contrast materi- 
als are subsequently employed to obtain finite order boundary/transition conditions. 
Finally, numerical results are presented in which the exact reflection coefficients are 
compared with those implied by the GIBC/GSTC to provide a measure of the con- 
ditions' accuracy and utility. 

In chapter three, the plane wave diffraction by a multilayer material slab recessed 
in a perfectly conducting ground plane is formulated and solved via the Generalized 
Scattering Matrix Formulation (GSMF) in conjunction with the dual integral equa- 
tion approach. This problem is significant in that a unique GIBC solution is obtained 
which can be used as a benchmark to test other GIBC/GSTC solutions. In the first 
part of the chapter we summarize the GSMF procedure. The dual integral equation 
method is then employed to formulate each of the subproblems and the neo"-sarv 
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solutions are obtained for both E~ and H z polarizations. These are given in terms of 
symbolic split functions which are then evaluated for the specific case of a multilayer 
grounded slab by casting the reflection coefficient in a form compatible with a GIBC 
simulation of chapter two. A number of scattering patterns are presented and the 

accuracy of the GIBC simulation is examined by comparison with known results for 
homogeneous slabs. 

In chapter four, the diffraction by a material discontinuity in a thick dielec- 
tiic/feirite slab is considered by modelling the slab as a distributed current sheet 
obeying generalized sheet transition conditions (GSTC). In the first section of the 
chapter, the GSTC representation of the distributed sheet discontinuity is used to 
develop dual integral equations in terms of the unknown spectral functions propor- 
tional to the sheet currents. These equations are then solved in the standard manner 
to yield expressions for the spectral functions in terms of unknown constants. The 
constants are dependent on the geometry and properties of the discontinuity, and 
are identified in this chapter for a few specific discontinuous layers whose diffraction 
solution is available. 

Chapter five deals specifically with the determination of the unknown constants 
for the solution presented in chapter four. This is accomplished by introducing a gen- 
et al eigenfunction expansion which is valid everywhere and subsequently recasting 
the solution obtained in chapter four into this format. This field is then analytically 
continued to the slab interior and continuity is applied at the material junction to 
provide the remaining constraints for determining the unknown constants. Specif- 
ically a point matching scheme is proposed in which an overdetermined system of 
equations is generated and solved for the constants using a least-squares technique. 
Various diffraction patterns are given validating the obtained solution for certain di- 



9 


electric half- planes 
are given for vario 


and metal-dielectric junctions on a ground plane. Finally curves 
us thick half-planes and junctions to illustrate their diffraction 


behavior as a function of thickness. 


CHAPTER II 


DERIVATION OF GENERALIZED 
TRANSITION/BOUNDARY CONDITIONS 
FOR PLANAR MULTIPLE LAYER 
STRUCTURES 


A GIBC/GSTC diffraction coefficient can only be as good or as versatile as its 
constituent GIBC/GS1C. Therefore, before diffraction coefficients of any generality 
can be developed it is necessary to construct GIBC/GSTC which are valid across all 
ranges of slab composition and thickness, and this is the task of this chapter. The two 
configurations considered herein are the multilayered slab having symmetric or non- 
s} mmetiic material composition (about its center) and the multilayered coating on a 
ground plane, as illustrated in Figures 2.1 and 2.2(a), respectively. The derivation 
of the GIBC/GSTC is accomplished via the Taylor series expansion method, whose 
versatility enables the treatment of non-planar as well as planar layers. In effect, 
the resulting conditions allow the simulation of the multilayered configuration as an 
opaque or transparent sheet (see Figures 2.2(b) and 2.3(b), respectively). 

In proceeding with the development of the GIBC/GSTC, we initially derive re- 
inrience lelations for the fields in adjacent layers. These are subsequently employed 
to develop infinite order boundary/ transition conditions that are conveniently ex- 
pressed as a matrix product. Approximations to the matrix element operators for 
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y 



Figure 2.1: Infinite multilayer slab. 
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2.3'. 'a) Zero thickness resistive and conductive sheet similation of multilayer 
slab, (b) Distributed resistive and conductive sheet simulation ot multi- 
layer slab. 
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lou and high contrast materials then lead to finite order boundary conditions. Fi- 
nally. numerical results are presented in which the exact reflection coefficients are 
compared with those implied by the GIBC/GSTC to provide a measure of the con- 
ditions accuracy and utility. 

2.1 Derivation of Infinite Order Conditions 


Consider the multiple layer slab with N v upper layers (y > 0) and \ L lower layers 
(U < 0) as illustrated in Figure 2.1. The physical parameters corresponding to the 
m ' n u PP er la - ver are denoted by a b m , & which refer to the relative permittivity 

and permeability, the index of refraction and the layer thickness, respectively. In 
a similar manner, the physical parameters corresponding to the m th lower layer are 
given by ^ *m’ T m- In the following we derive transition conditions to effectively 
replace their presence with a distributed current sheet. We begin this derivation by 
first introducing a relation between the fields on the two sides of a single layer. This 
is generalized to relate the fields of distant layers and those at the upper and lower 
boundaries of the fictitious current sheet. For convenience, we may consider both 
polarizations simultaneously by introducing the definitions 

E y , E y polarization ( H y = 0) 


F, = 


u = 


Hy, Hy polarization (E y = 0) 
e, E y polarization 


12. T 


( l-i, Hy polarization 
l sing a Taylor series expansion, the normal field components at the top and 
bottom of the m th upper layer may be related as 


10 1 ikr V 

“mi r y J y=7 c + r V + ... +(r V_ i)+ = U m E — y=rU + ... + T U_ i+{r u ) - 




X p I _ V' -l + 1 r . 

^ y- r i +■■ •+( r m_ 1 ) + 2 -* d y *yjy=Tt ’+ T 2 +-+T%-i+(r&)-- '--dj 



15 


where 


T ± = T ± 0 + 


(2.4) 


and F is an operator defined as 


- J d 
by ~ kdy 


(2.5) 


Clearly. (2.2) - (2.3) provide relations of the fields within a single layer and if we are 
to derive a condition relating the fields at the top and/or lower surface of the slab, 
it is necessary to establish similar relations among the fields in different layers. As a 
first step towards this goal we may proceed to express the right hand sides of (2.2) 
and (2.3) in terms of the fields above the boundary y = . To do so. 

it is instructive to resolve the right hand sides of (2.2) and (2.3) into a summation 
of odd and even derivatives of F y . Subsequently, the wave equation may be invoked 
to rewrite the normal derivatives in tangential form, thus allowing the application of 
the field continuity conditions. We have 


*V = 


*1 


( 2 . 6 ) 


with 


% = K + K 


w 


here 


8* = f-- 


6? = 


Li- 
fe dx 

Li- 

fe dz 


(2.S) 


Invoking now the continuity of u m F y and & l y F y (including their tangential derivatives) 
across the layer boundaries, (2.2) and (2.3) may be rewritten as 


F\y=r e + ,e + ... + (r v_ l) + = Qu(u L L K H^f) • F\ y=T u +T u^ Mr u^ 


(2.9) 
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In this 


F = 


Q C (U. K.T.Sf) = 


uF j 

S y r j 

qil(u, K , T. 8 {) q 12 { U. K. T. <5 £ 2 ) ' 

92l(«, «, 7-, S‘) q 2 2 \U, K, T, 81) „ 


2 . 10 ) 


( 2 . 11 ) 


^n(i/,«,r,(5 ( 2 ) = cos (k 0 rjttj 

k oT ^rr^ 


sm 


( h 2(u. K.r, = ju- 


\/' t2 - *? 


V 2 i( w, k, r, <5;) = - ^7 sin ^’ 0 r^ 2 - £j? 

<?22( «, T, Si) = cos (^k Q T\flF - £ f 2 j . 

Expressions relating the normal fields in adjacent lower lavers mav bt 
in a similar manner. We have. 


f\ 


^ = ~ r r~ r 2 


r m-l) + 




ri-r/-..-(r4)+ 


( 2 . 12 ) 

obtained 

(2.13) 


Q l {u,k.t.81) = 


(2.14) 


where now 

q n (u,K,TJl) -qi2{u,K,T,Sl) 

\ -<? 2 l(u, K,T, 8f) q 2 2{u,K,T,8'1) J 

Equations (2.9) - (2.14) constitute fundamental recurrence relations for devel- 
oping multilayer GIBC/GSTC. Each of the <711,912, etc. is an infinite-order linear 
differential operator in even powers of <5 2 . This is evident when the sin and cos terms 
are cast in their Taylor series representation (note that the square root functions ap- 
pealing in (2.12) do not have branch cuts). The finite-order boundarv conditions are 
then derived by truncating the Taylor series representation of 911,912,921, and q 2 2 - 
Applying (2.9) and (2.13) recursively, we may establish a relation between the 
fields at the top and bottom of the layer. We have 


' Ah 

n<? 

m = 1 


K m- T m 


L X;) 


' F}y=-T 1 L -Tf--,..-{T* i )+ — -^Jy=0- 


. 15 ) 
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.w 

no 

•n = 1 


C‘( U !n- K m > ' / 


■ F] y= .u^ r c + Mr u )+ = F j,=c 


. 2 . 1(1 


and since F is continuous across y — 0. (2.15) - (2.16) imply the GSTC condit 


ion 


.v 


L _ 


I ] 


m = 1 

Xu 


' F\ y =_ T L_ T l_.„_ lT L 


( 2.1 


n Qu( U m' 


771 = 1 


* y = r 1 L ' r + r^ + ... + (r^) + 


-Cl 


The part of the slab occupying y u < y < rf + r? + . . . + r^ r may now be replace 
bv the upper background medium (with its geometrical and material parameter 
denoted bv the subscript b and the superscript U) while the part of the slab occupying 
y L > y > — — rk — . . . — r y is replaced by the lower background material (with 

geometrical and material parameters denoted by the subscript 6 and the superscript 
L). Using a Taylor series expansion the boundary fields may be related to the 
equivalent fields at y = y u and y = y L . In so doing, we obtain 




771 = 1 

Xu 


Qu( U b i K b i T b ’ -U.jJ y=y L 




771 = 1 


Q l( U b > K b i~b i ) F e,J y=y u • 


(2. IS) 


where 


Nl 


= Y.rt-{-y L ) 


i=i 

Nu 


r? = EF -( y u ) 

i=i 


(2.19) 


Although compact, the transition conditions (2.18) provide little insight into the 
phvsics thev represent. It is therefore instructive to reorganize them in a form that 
leads to its physical interpretation. To this end we introduce the definitions 

r ^ 2 ) £; 2 (*?n 


.Vr 


n q 


, r 7l= 1 


n < -rL f2\ 

Qu( e b » K b > T b ) — 


CM) C% 2 {6\ 


t ) > 


-H + +H ■’?-= -H + , 



and i\F^ f are given as 
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in which 


St 

def 

5 r x + 6 z z 

F ; 

riof 

E r x + E y y 

H t 

def 

H x x + H y y 

E t 

def 

E? W 

E y 

def 

£,*'W 

u; 

def 

HT W 

n; 

def 

H? W- 


We also note that the superscript notation e and h refer to E y and H y polarized 
excitations, respectively. 

In view of (2.24). the transition conditions given by (2.21) and (2.22) are now 
readily interpreted as a representation of "distributed” resistive and conductive sheets 
occupying the volume y u < y < y L and supporting equivalent electric and magnetic 
currents (see Figure 2.3a). These, of course, give rise to discontinuities in the fields 
(and their tangential derivatives) at the upper and lower boundary of the sheets. The 
equivalent electric and magnetic currents are denoted by J e ' h and M e,h , respectively, 
and are defined in (2.24) above. Additionally, the subscripts J and M appearing 
in (2.23) denote field discontinuities traditionally associated with the presence of 
electric or magnetic currents, respectively. For convenience, the currents and result- 
ing discontinuities A F are presented in terms of both normal and tangential fields. 
When y l = 0 + and y L = 0“ (see Figure 2.3b), the distributed resistive and con- 
ductive sheets are "compressed” onto an infinitely thin sheet occupying the plane 
y = 0. Such thin sheet representations are attractive for the application of transform 
techniques in diffraction problems and are generalizations of the resistive-conductive 
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sheet simulations given in [59], [39] and [31]. In contrast to the referenced simu- 
lations. the ones given here exhibit coupled electric and magnetic currents, whose 
presence complicates the application of analytical techniques. This coupling is at- 
tributable to the distributed nature of the polarization currents across the original 
slab configuration. If, however, the multilayered slab is symmetric about y = 0 and 

U = ~ V L • then it can be shown that M and J decouple. In particular, for this 
special case 


C\\ — 

£u = 

C 2 1 = 
C-22 ~ 

leading to 

u’ u (S?){af;} = 

= 

UfoSfHAFh} = 
U^SlHAF}} = 


till 

—U 12 
—U 2 i 

^ 22 , (2.26) 

-w l2 (8 2 t ){r} 

(2.27) 


when substituted into (2.21) and (2.22). Evidently, the coupling of the current 
components depends on the degree of asymmetry in the slab and an assessment on 
the level of coupling can be obtained by comparing the magnitude of the operator 
coefficients in £n — U n , £ 12 -fi Ui 2 , £21 + £Ci and C 22 ~ W 22 relative to those in 
Oi + r’Y u , £ 12 — U\i, £21 — U 2 \ and £ 22 + U 2 2 . 

If a ground plane is inserted in the symmetric slab at y = 0, the resulting structure 
becomes opaque and its sheet simulation is further simplified. For the case where the 
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ground plane is a perfect electric conductor (PEC), the electric currents are shorted 
out and the second and third equations of (2.27) become the boundary conditions 
corresponding to a coated PEC. Conversely, a perfect magnetic conductor (PMC) 
shorts out the magnetic currents, leaving the first and fourth equations of (2.27) as 
the corresponding sheet simulation. An alternative approach to deriving the GIB(_ 
corresponding to coatings is to employ image theory in conjunction with (2.27). The 
field components in the equivalent image configuration of a PEC grounded slab are 
related as 


E y {y> 0) = E y {y < 0) 

6yEy(lJ > 0 ) = -8yEy{rj < 0 ) 

Hy{y > 0 ) = -Hy{y < 0 ) 

8 y H y (y> 0 ) = 6 y Hy(y < 0 ) (2.2$) 

and when these are subsequently introduced in (2.27) vve obtain the boundary con- 
ditions 

+ uM = °- <-*» 

corresponding to a slab on a PEC. Similarly, the image fields for a slab on a PMC 
satisfy the relations 

E y (y > 0) = -E y (y < 0) 

8 y Ey(y > 0 ) = 8yEy(y < 0 ) 

H y (y> 0) = H y (y< 0) 

8yHy(y > 0 ) = -8 y Hy(lJ < 0 ) 


i 2. 'll) i 


oo 


which lead to 


$ir x ^& 2 t ){E+) +ir n (5 2 t ){6 y E+) = 0 

tfU}xtf){Ht}+Uk(6t){6 y H+} = 0. (2.31) 

when substituted in (2.27). Expressions (2.29) and (2.31) represent opaque conduc- 
tive and resistive sheets, respectively, and are the dual of each other. We note that 
it A t - = 1. (2.29) reduce to those given in [3 1 J . 

' •' - uinaiize the above development, sheet simulations were derived that model 
(or replace the presence of) the multilayer slab and coating. Equations (2.21) through 
(2.27) are referred to as generalized sheet transition conditions (GSTC) for the 
transparent resistive/conductive sheet representation of the multilayered slab. On 
the other hand, equations (2.29) and (2.31) are described as generalized impedance 
boundary conditions (GIBC) f-" the opaque sheet representation of the coated ground 
plane. They aie given in a compact matrix form and are valid for anv arbitrarv finite 
number of layers. Their versatility, however, is offset by the presence of infinite order 
derivatives as implied by the definition of the operators, thus, limiting their applica- 
bility to analytical and numerical treatments. It is therefore appropriate to consider 
finite order approximations of the operators leading to conditions of practical use. 
In the following we consider such approximations of the operators on the assumption 
of low contrast (small k) and high contrast (large k ) layers or coatings. 

2.2 Low and High Contrast Approximations for Matrix El- 
ement Operators 

Low contrast approximations to the matrix element operators may be derived 
by leplacing the trigonometric functions in (2.12) by their Taylor series expansions 
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which are then truncated. To 0(r u ) these become 


/here 



int{¥) 

m!(f) 


i]ll \ U. K . T. df) = fj22{ 11 • K - < 

-7)= E 


k , 7. /. m) 


(=0 

m = l 





(«?)' 

,/ i: ( u. k. r. 6 2 ) % ju Y 

Y a 2 {K,T,Lm) 

1=0 

j m = l 



■ int[X± 1) 

mK^) 

~ 

/ V / 

q 2 l(u. K.T, 8?) « - Y 

Y 03(K 

, r, /, m) 

(^ 2 ) m 

u /=o 


- 


n .( ur T 1 rn\ — 

l) ,+m (m)l(A: 

_ ■ j 2 m 2 ( m - 1 ) 




r, /, m) = 
a3(«,r, /,m) = 


(m-/)!(/)!(2m)! 

■l)' +m (rn)!(fc 0 r) 2m+l /c 2(m -' ) 


(m - /)!(/)!( 2 m + 1)! 
:-l) ,+m (m)!(fc 0 T 




(m - /)!(0!(2m - 1)! 
It is a simple matter to demonstrate that the substitutions 


\ji^ 2 - «6 = f * 




f 2.3: 


[2.33 ] 


(2.34) 


applied to (2.32) lead to a normal derivative representation of q n - qn- 

To obtain high contrast (large k) approximations for q t] , it is necessary to utilize 
the binomial expansion 

M) 


yjt 2 - Sf = K 


1 _ I (hY - J_ (h. V - L ~ 3 

2 \ k / 2 • 4 V K ) 2-4-6 J 


< 1. 

(2.35) 


Substituting this into (2.12) and again employing a Taylor series repiesentation tor 
the trigonometric tunctions, we obtain to 0(k '^) the appioximations 


r/n(u, k , t , 6 2 ) =2 cos(A.- 0 r«) 
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a; 

-E 

p=i 


A/ 

E 

/=max{ 1 ,2p — Af } 


r 1 (U.r)T(/.O.p) 

K'P-‘ 


(«?)' 


>( u. K, r. <5;) 


— sinf k 0 r k) 

K 


Af-l 

+j u E 


A/ - 1 

£ 


P=1 L^=77iar{1.2p-A/-f-l} 


r 2 (/, /c. r )T ( /, 0. p) 
k 2 p+ 1 - 1 


\r-2 

-ju E 


A/ — 2 

£ 


'Z’lO- K. r, 6j 


where 


p— l [_/=mar{l ,2p— A/-h2} 

= — sinlA,’ 0 r«) 

u 


r 2 (/ - 1, /t, r)T (7 - 1. Lp) 
^2^+2-; 


j -V+i 

-i V 

, 2-w 


E 


r 2 (/.K.r)T(/,0,p) 


t>=l [_/=max{l ,2p — A/ — 1 } 


K 2p-l-l 




--£ 
U , 
p = l 


A/ 

E 

/=m ax{ 1 ,2p— A/ } 


r 2 (/- l,K,r)T(/,0,p) 

K 2p-l 


-7\P 


(«?)' 


(2.36) 


T x (I.K,t) 


{-iy-t(i/2) { k oT y 

' COS (k 0 TK) 


(/)! 

k sin (Ic 0 tk) 

T 2 ( /, K. T ) 


(_l)<nf<// 2 )( fcoT )/ 

f si n(k 0 TK) 



(/)! 

k — COS ( k 0 TK 

T (l x J 2 .p) 

= 

E bi{ii)bi{i 2 ) • • • b x {ii 1 )b2(ii 1+x 

b i{p) 

= 

(2 p - 3)!! 

2 p p! 


hip) 

= 

(2p - 1)!! 
2 p p! 


X\\ 

def 

X(X-2){X-4) 

■ ■ ■ (3)(1) ; 


/ is even 


tie/ 


(2.37 


The sum defining the function T (/j,/ 2 ,p) includes all product terms satisfying the 
relation i m = p. It should also be noted that l x and / 2 denote the number of b x 
and b 2 . respectively, comprising the product terms. For example, if l x = 2, l 2 = 1 
and p = 3 then T(2,l,3) = 6 1 (1)6 1 (1)6 2 (1); if l x = 2, l 2 = 1 and p = 4 then 
T ( -• 1 • 4) = b x (2)b x ( 1)6 2 ( 1) + 6 1 (1)6 1 (2)6 2 (1) + ^i(l )6 1 (1)6 2 (2) and if l x = 2, /? = 0 
and p = 4 then T (2.0.4) = 6 1 (1)6 1 (3) + 6 1 (2)6 1 (2) + 6 1 (3)6 l (l). Corresponding 
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Order 

Low contrast q tJ coefficients 

_o 

7 1 

'7 11 = q:2 = 1- <7 12 = <721 = 0 

_1 

' 1 

<7ll = f{22 = 1* 4 ll'2 ~ J u \k r \' c l2\ = ^ T l 

2 

7 f 

<711 - <722 - 1 - ■^t-(k 2 <$f ) • 712 = juikr u q 2 \ = ^t x {k 2 - S') 

_ 3 

' 1 

<711 - 722 - 1 - ( \ 2 1 ' (k‘ 6;), 7l2 - jUlkTi [l (k‘ G 2 )]' 

, s ,=jtr,( K ’ -<?)/«, [l-^(<c 2 -«?)] 

' 1 

= qn = 1 - ( K 2 - 6*) + («* - f!f 

721 - jbT X (K- - 6;)/u x [l-^(x 2 - 6f)} 


Table 2.1: Low constrast approximations to q tJ operators. 

expressions in terms of normal derivatives can also be obtained by employing (2.34) 
in (2.36). 

Expressions (2.32) - (2.37) represent finite order approximations to the operators 
<7n- <7i-2,92i! an d <722 1 and can be used to generate finite order GIBC or GSTC. To 
do so for a specific multilayered slab or coating, one first examines each constituent 
laver and approximates the corresponding matrix elements by their low contrast 
expressions (2.32) - (2.33) or high contrast expressions (2.36) - (2.37), as appropriate. 
As an example, basic low contrast and high constrast approximations to the q tJ 
operators are give in Tables 2.1 and 2.2. These finite order expressions are then 
substituted into the matrices Q v and Q L given by (2.11) - (2.14). The resulting 
expressions for the IA and C matrix operators are finite polynomials in even powers 
of 6‘. In passing, we note that the simulations presented here-in may also be extended 
to model a longitudinally inhomogeneous slab once this is approximated as a layered 
slab such as that shown in Figure 2.1. 
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Order 

1 

1 High constrast < 7 , ; coefficients 

«° 

9n = 922 = cos (kIctx). q u - 0 .< 7 21 - ^ sip/.vfrrj - C os(/cfcr x ) 

K~ l 

( hi — 922 cos (/cti t ) + sin ( k k . q v2 — ^ sin (Kkr x ) 
<721 u * siii(kA:/x) {l j{r ■ ' } -j^-cosf/cArr!) 

hi “ 

9u 922 — cos (kIcti) jl — 1 | + sin (kIcti ) — 

( 7i2 ^ sin (hckii) — ~ 2 1 cos (kIcti ), 

K-- 3 

9n = 922 = COS {nkTi) {l - } 

+ «„(*tn){ij4 + iSt-^5fi}, 

In = sin (Kkn) {l + £ - lijjja} - cos (KkTi)t 

921 - i £ sin(ic*r 1 )|l & + 1*3 J 4 *? \ 

1 “1 V iJ l 2* 2 8* 2 8k* 1 388k* )■ 

-^cos( Kjfcn)/^ *il£ 

A/ L 2* 8x 3 48« 3 /’ 


Table 2.2: High constrast approximations to q tJ operators. 
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Figure 2.4: Coordinate system for reflection and transmission coefficient derivation. 

2.3 Evaluation of the Boundary/Transition Conditions 


To evaluate the accuracy of the derived boundary/transition conditions one ap- 
proach is to compare the plane wave reflection and transmission coefficients implied 
by the finite order sheet simulation with the corresponding exact coefficients. 
Consider the plane wave 




Hi 


E. 


yo 


0 ]kK% ( x cos a sin/ 3 +y sin a +2 cos a cos 0 ) 


(2.38) 


H 


yo 


incident upon the sheet satisfying a given GIBC or GSTC (see Figure 2.4). The 
generated reflected field can then be written as 

O 

ReE v q 


F T 
^ y 


HI 


x cos ot sin 0 — y sin a + 2 cos ar cos 0 ) 


RrH 


H™ yO 


(2.30) 


2S 


and in the case of a GSTC, the transmitted field takes the form 

f Te EyO 


E l 


^ ( x cos sin h?' + y sin oe f z cos ot f cos 3 * ) 




1'h HyO 


(2.40) 


To find the reflection coefficients Re,h and corresponding transmission coefficients 
1 E.H ue substitute (2.3b) - (2.40) into one of the boundary or transition conditions 
given by (2.21). (2.22), (2.27), (2.29) or (2.31). By carrying out the differentiations 
in a straightforward manner we find 

multilayer slab 
multilayer symmetric slab 


Re = 


*?m« ~ R l'c 


i-R e + 

2 pec r 2 pmc 


R 


pec 


T e = 


De 

^pmc 
^pec ” ^pmc 

if? e _ IP* 

2 A pec 2 pmc 


multilayer coating on PEC 
multilayer coating on PMC 

multilayer slab 
multilayer symmetric slab 


multilayer slab 

^*-pec ^pmc 


Rh = l 


T h = 


In (2.41) - (2.44), 


1 t 1 ph 

2 l pmc * 2 ^pec 

R h 

C pec 

r>K 

^pmc 
R pmc ~ R pec 

- R h — 4 R h 

2 1 l pmc 2 " rl pec 


multilayer symmetric slab 
multilayer coating on PEC 
multilayer coating on PMC 

multilayer slab 
multilayer symmetric slab. 


R e = 

pec 


^ pmc 


R e - 


(2.41) 


(2.42) 


(2.43) 


(2.44) 


k[ t sin a U\ 2 

(k^ ) 2 cos 2 a 

- 

(«^) 2 cos 2 O’ 

sin a 

sin q U 22 [{ *■¥ ) 2 cos2 a . 
sin a li\ 2 f(«f ) 2 cos 2 a] 

I + e 6 ^21 1 
- 

(4) 2 cos 2 a] 

[(Kj') 2 COS 2 Q 

sin a 

k ^ sin a U{ 2 
p sin a' C% 2 

[(4' ) 2 cos 2 a] 
(k \ ) 2 cos 2 a 

+ 4'^u [ 
] - e 6^21 

(4 ) 2 cos 2 a] 
(«f ) 2 cos 2 a 

. „jk sina'4-^f 

k[‘ sin q U 22 

(4) 2 cos 2 a 

+ ( 

(k^) 2 cos 2 a 

— g L 0 ^ 0 


■yV' 
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4 sin a' C\ 2 («5 ) 2 cos 2 q 

r"> £ 1 - J 

- t[C e u [(44 cos 2 a 

y L sin ot ' sin a] 

pmc “ K l h sill a U e n [(4) 2 cos 2 a; 
4 sin a lih l ( 4 ) 2 cos 2 a 

ph ' 1 

1 + 4 U u [(44 cos2 a ] 

- 4 b U\\ [<'44 cos 2 0 ] 

. e J- k «b'y u sina 

V ' ; ~ 4 <ln a l4 : [( 4 )- cos 2 a\ 

4 sin a [(4 ) 2 cos 2 0 ] 

V 4 4l [(.4 r COS" Q ] 

- 44h [(44 cos 2 a] 

_^j2kK^ y u sin a 

)pmc sin a U -22 [(*& ) 2 cos2 a \ + t- l b^ 2 \ [( K i )“ cos " Q ] 

= „ 4 sin a' £* 2 j(4’) 2 cos 2 a] - 44i [(44 cos 2 a 

l e J k [*by L sinck' + .K^ y 1 ' sin a] 

pec af sin a 14 [( 4') 2 cos 2 a 

4 sin a' £ 22 [(4 ) 2 cos2 a 

ph *“ l 

] + 4^l\ [(44 cos2 a ] 

] - 4^21 [(44 COS 2 Q 

]_.?«[«£ V L sinar'+.K^y 1 ' sine* 

pmc ~ 4 sin a 14 [( 44 cos 2 a 

i + 4'Un [(44 c °s 2 « 


(2.45) 

an from Snell's law 

4 sin a = y/(4) 2 - (4) 2cos2a - (2.46) 

To obtain a composite sheet simulation of a multilayer slab or coating, it is 
necessary to first model each layer individually through its Q v matrix. The matrices 
are subsequently combined to yield the £ and U operators of the composite boundary- 
condition. Thus, the accuracy' of the overall simulation can be assessed by examining 
that of the individual layers comprising the slab or coating. In the case of a single 
layer simulation (with non-shifted surface) the £ and U operators reduce to £n 
£ 22 = 1, £ n = £21 = 0 and U, } = q, r This simplifies the analysis and in the 
following we examine the accuracy of the proposed GSTC simulation as a function 
of the condition’s order. Only the £ y -polarization incidence is discussed but similar 
results apply to the //^-polarization case as well. 

Figures 2.5, 2.6, 2.7 and 2.8 present the maximum error in \R ex — R gs tc\ over 
real angles as the layer thickness is varied. The data in Figure 2.5 corresponds to 
a low contrast simulation with e r = 2 and fi r = 1-2 as the order of the transition 
condition increases from 2 to 9. In the region where |/7 ex Rgstc\ ^ lcSS tll&Il .15 
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Figure 2.8: E y polarization high contrast simulation of a single layer having 
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creases. This is indeed observed in Figure 2.S where the layer's relative constitutive 
parameters are e r = 5 and fi r = 3. As shown, a layer of up to 2 A thick can be 
accurately simulated with a 12 e/l order condition. Figures 2.5, 2.6, 2.7 and 2.8 are 
recomputed for H y polarization in Figures 2.9, 2.10, 2.11 and 2.12. respectively. The 
results are very similar to the E y polarization case, and the same remarks made 
above are applicable here. 

Using the single layer data, such as those presented in Figure 2.5-2.12. it is pos- 
sible to synthesize a. multilayer simulation. As an example, consider an E y polarized 
plane wave incident on a three layer slab having = 5 — jO, = 3 — jO, = .4 
for the first layer, = 3 5 ~ j0, = 2 - j 0, rf = .4 for the second layer and 

= 2 — j 0, P 3 = 1.2 — jO, 7$ — .2 for the last top layer (see Figure 2.13). To 
select the individual layer models for E y polarization we conjecture that the maxi- 
mum error of | R ex — R g 3 tc\ for the composite sheet will be bounded by the sum oi 
the maximum errors of \R ex — of the constituent layers in isolation. Examining 

Figures 2. 5-2. 8, we observe that the sum of the errors of the 6 th order high con- 
trast representations of the bottom and middle layers and the 5 th order low contrast 
representation of the top layer amounts to .062 as required by the design criteria. 
This suggests that the composite sheet employing these representations will have an 
acceptable performance when the total thickness is 1A. The actual maximum error 
of the designed simulation is compared with the sum of the isolated layer errors in 
Figures 2.13. The corresponding errors with a PEC and a PMC inserted at y = 0 
are also plotted in Figures 2.14 and 2.15. We note that the sum of the isolated layer 
errors serves as a reasonable upper bound for the new simulation design. 
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Figure 2.9: H y polarization low contrast simulation of a single layer having relative 
constitutive parameters e = 2, u = 1.2 
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2.4 Concluding Remarks 

Tliis chapter has dealt with the development of generalized impedance boundary 
and transition conditions for multilayer coatings and layers. A major effort was also 
devoted to present them in a form which is compact and convenient for further usage. 
Because of the generality of the derived conditions they are ideally suited lor use in 
the following cl voters where general diffraction coefficients are developed. 


CHAPTER III 


DIFFRACTION BY A MULTILAYER SLAB 
RECESSED IN A GROUND PLANE VIA 
GENERALIZED IMPEDANCE BOUNDARY 

CONDITIONS 


The canonical geometry formed by two semi-infinite planar slabs joined end-to- 
end is a difficult one to model properly, particularly if the slabs are thick. For this 
reason very few computed results exist to verify GIBC/GSTC diffraction solutions 
obtained for thick slab junctions. This need is addressed herein by computing a 
unique alteinati\e GIBC/GSTC solution for a specialized thick slab geometry. This 
solution provides a benchmark which can be used to partially verify more general 
GIBC/GSTC diffraction solutions. The special problem considered in this chapter 
is the plane wave diffraction by a multilayer slab recessed in a ground plane in 
Figure 3.1, and the alternative method used is the generalized scattering matrix 
formulation (GSMF) [25]. 

Related but simpler geometries have been studied in the past [2] [6] [20] [11] [26] 
[*^] 1^1 [— b], and most of these solutions involved the interior and exterior 

fields. This is alleviated herein by modeling the slab as a surface characterized bv a 
plane wave reflection coefficient R ( cos </> 0 ) (see Figure 3.1), where <f> 0 can be extended 
thiuugn anal) tic continuation in the complex plane. This enables us to carry out 
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(b) 

Figure 3.1: (a) Multilayer slab recessed in a PEC ground plane, (b) Representation 
of slab as surface with reflection coefficient R. 

the analysis in a symbolic manner regardless of the inhomogeneity profile of the 
layer. Hence, although our focus in this chapter is the multilayer recessed slab, the 
derivations will be applicable to any vertically inhomogeneous slab. 

The problem herein is formulated via the dual integral equation approach [10] in 
conjunction with the GSMF. The GSMF is applied to the recessed stub structure 
of Figure 3.2(a), depicting a perfectly conducting half plane elevated a distance 6 
above a reflecting surface with a perfectly conducting stub recessed a distance cl 
away from the half plane edge. This formulation requires the solution to a number of 
individual subproblems. As illustrated in Figure 3.2(b)-(f), they correspond to the 
problems of direct diffraction, mode coupling, mode reflection, and mode radiation. 
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(f) 


Figure 3.2: Illustration of recessed stub geometry (a) and associated subproblems: 
(b) direct diffraction, (c) mode coupling, (d) stub reflection, (e) mode 
reflection at the waveguide mouth, (f) mode launching. 
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Once solutions to each of these subproblems have been obtained via the dual integral 
equation method, they can be combined in accordance with the GSMF prescription 
to yield the diffraction for the original structure in Figure 3.1(a). 

Unfortunately, the conversion of the symbolic solution into one of practical use 
proves to be a formidable task when the reflection coefficient of the grounded slab is 
obtained in its exact form. The fundamental difficulties are related to: 

• the factorization or splitting oi the associated Weiner-Hopf functions into com- 
ponents regular In the upper and lower half complex plane, and 

• the extraction of the complex zeros (i.e., the waveguide modes) associated with 
the split functions. 

The pertinent Weiner-Hopf functions cannot be factored analytically and one must 
therefore resort to a numerical scheme (e.g., see [28]). Also, in solving for the complex 
roots of the pertinent split functions, it is necessary to employ a search algorithm in 
the complex plane, a process which is numerically intensive. We circumvent these 
difficulties by replacing the grounded slab by an opaque sheet satisfying a GIBC of the 
form given in chapter two. Under a GIBC approximation, the approximate reflection 
coefficient is cast as a ratio of polynomials in cos (f> or sin making the determination 
of the complex poles and zeros of the reflection coefficient a simple task. As a 
result, the required W'einer-Hopf factorizations can be obtained analytically leading 
to computationally efficient solutions. 

In the first part of the chapter the GSMF procedure is summarized. The dual 
integral equation method is subsequently employed to formulate each of the sub- 
problems and the necessary solutions are obtained for both E z and H, polarizations. 
These are given in terms of symbolic split functions which are then evaluated tor the 
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specific case of a multilayer grounded slab by casting the reflection coefficient in a 
foini compatible with a GIBC simulation ot the slab given in chapter two. Results 
are given and the accuracy of the GIBC simulation is examined by comparison with 
known results for homogeneous slabs. 

3.1 Description of GSMF Procedure 


In this section, the generalized scattering matrix formulation (GSMF) is applied 
to the geometry given in Figure 3 ‘1(a). This consists of a perfectly conducting half- 
plane located a dist? ce S above the grounded slab, with a perfectly conducting stub 
recessed a distance d away from the half-plane edge. To concurrently treat both the 
E : and H . polarizations ot incidence, the quantities F. and F x are introduced. They 
are defined as 


F, = 
F r = 

and from Maxwell’s equations 


E„ 

E z polarization, 

Z 0 H z , 

H z polarization. 

Z 0 H X , 

E z polarization, 

E x , 

H z polarization. 


(3.1) 

(3.2) 


F x = 


jjh dF z 
k dy 


(3-3) 


where 


Uj = 


(3.4) 


— 1, E 2 polarization, 

l 1, H z polarization. 

The individual problems to be considered in the GSMF prescription [25] are as 

follows: 


1. Evaluation of the direct diffracted field by the substructure in Figure 3.2(b) 
due to a plane wave incidence. This field can be expressed as 

F? d (<f>, <Po) = f Pdd (cos a, cos 0 O ] 6) e~ jk ° lJCOsl ' a ~' t>) da TV 

J c 
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where P dd is the spectrum associated with the currents induced on the half 
plane and {p. o) are the usual cylindrical coordinates of the observation point. 
Additionally. C' is the complex contour composed of the directed line segments 
:0 — j dc. 0 — j Oh [0 — j 0, - — jO], [x — j 0, x + joe] in the complex a plane. 

2. Evaluation of the field coupled into the waveguide due to a plane wave incidence 
(Figure 3.2(c)) as illustrated. We denote the field corresponding to the n th 
coupled mode as 

F c ;n (• <p 0 ) = C n (cos 0 O \ S) e~ jknX (3.6) 

where C n (cos 0 0 \ 6) is usually referred to as the coupling coefficient and k n is 
the propagation constant associated with the n tfl mode. 

3. Evaluation of the modal field reflected at the stub (Figure 3.2(d)). This can 
be expressed as r mn e jfcmr where T mn is the stub reflection coefficient of the n th 
mode to the m th mode. 

4. Evaluation of the reflected field at the waveguide mouth due to the n th mode 
(Figure 3.2(e). This can be expressed by R mn (<!>) e~ jkmX where R mn (<5) is the 
reflection coefficient of the n th mode to the m th mode. 

5. Evaluation of the radiated field attributed to the m th mode incident at the 
waveguide mouth (Figure 3.2(f). This field can be expressed as 

F l zm {0)= f P m (cosa-6)e-^ cos ^ ) da (3.7) 

J c 

where P m (cos a, if) is proportional to the spectrum of the currents induced on 
the half plane due to the incident m th mode. 
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Accordingly, the scattered field by the recessed stub geometry in Figure 3.2(a) is 
given by (for ij > 8) 

F’(o.o 0 ,6.d)= / [P dd (cos a, cos o 0 \ 8)+P mod (cos a, cos d 0 [ 8. d)l e ~ jkofiCO& (a " d b/a 
J c 

(3.S) 

u here P m od (cos a, cos <p 0 ] 8 , d) is associated with the presence of the stub and includes 
the contribution of the waveguide modal fields. It can be written in a matrix form 
as [25] 

Pmod (cos a. cos o 0 : 8, d) = 

[■Pm (cos a; £)] 7 {[/] - (rf)][r TO „][PV TO „ 

'[^mn (d)][r mn ][IF mn (d)][C n (cos do] d)j (3.9) . 

in which the brackets signify column or square matrices depending on whether one 
or two subscripts appear, respectively. In addition, [/] denotes the identity matrix 
and [H mn (d)] is the modal propagation matrix whose elements are given by 

f e~ jknd , m = n 

II mTj (d) = < (3.10) 

l 0, m ^ n. 

To obtain the field scattered by the recessed material slab it is only required to set d 
and 6 to 0 in (3.8) and (3.9). In this case, [ W mn ( d = 0)] becomes the identity matrix 
and [T mn ] reduces to [/] or — [/] for H z and E z polarizations, respectively. Thus, 
Pmod becomes 

Pmod (cos ct, cos <Po) *= P mo d (cos a, cos do] 8 = 0, d = 0) 

= [Pm (cos a )]{[/] — Vl[/? m n]} 1 [v\C n (cOS <po)) (3.11) 

where 


def 


P m (cos a) 


P m (cos a; 8 = 0) 
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C n (cos Oq) d - c n (cos <? c ; <5 = 0) 

Rmn = Rmn(S = 0) . (3.12) 

Hence, the field scattered by the vertically inhomogeneous recessed slab can be ex- 
pressed as 

F*(o.6o) = /P ra (cosa.coso 0 )e-— ; y > 0 (3.13) 

where 

P ra (cos a. cos o 0 ) = (cos a. cos 0j) + P mo i (cos a. cos o 0 ) . (3.14) 

The steepest descent method can then be employed to evaluate (3.13) and obtain 
the diffracted field. 

3.2 Plane Wave Diffraction and Mode Coupling 


Consider the plane wave 

pi __ e jk 0 (xcos<t > 0 + y sinrio) 

F' x = u 1 sino 0 e jMxc “*° +ysin ' M 


( 3 . 15 ) 

( 3 . 16 ) 


incident at an angle <p 0 upon the structure depicted in Figure 3.2(b). In the absence 
of the perfectly conducting half-plane, the total fields may be written as (for y > 0) 


ppvu 

= fi + f: 

( 3 . 17 ) 

Fr 

= f' x + f; 

( 3 . 18 ) 


with 


FI = R{cos<p 0 )e jk ° (xcosrt> °- ysin * o) (3.19) 

F r r = -u 1 P(cos<p 0 )sin<p 0 e J ^ (xcos ^- ysm ^ ) (3.20) 


where R( coso 0 ) denotes the plane wave reflection coefficient of the grounded slab 
telened to y — 0. For the general case of a plane wave incident at an angle a on a 
\eiticallv inhomogeneous ground plane coating, R ( A = cos a) may be represented as 


R(\) = - [ - 4 ( A )-Vl-A 2 g(A) ' 

LA (A) + ' 

where -4(A) and B (A) are even functions of A, with any branch cuts in 
appearing in B (A) and vice versa. 

The introduction of the perfectly conducting half plane at y = 8 ge 
additional scattered field component F’ so that the total fields become 


(3.21) 
.4(A) also 


nera* an 


F z = Fr + F: (3.22) 

F x = Fr + F’ (3.23) 


This scattered field is due to induced currents on the perfectly conducting half-plane, 
and can thus be represented by an angular spectrum of plane waves. A suitable 
representation is [ 10 ] 



(f c P (cos a) e jko($ sin a e~~ jkop cos da 


y > S, 


f: = < 

f c Q (cos a) e-i koSsina | e “^^ C08 (^+ 0f ) 



( 3 . 24 ) 


+ R (cos ( 7 T — a)) e -jkoPcos((t)-at) 

da 

0 < y < 8, 


implying 






f - fc v 1 sina P (cos a) e ]koSsina e~ jk <=/ ,<:os ^-o) ( la 

y > 6, 


ii 

r 

fc 1 ? 1 sin a Q (cos a) e jkoSsinot jk a p cos (<t>+ot) 



(3.25) 

i 

— R (cos (tt — q) ) e _Jf;Cop cos( 0 -Qf)J 

da 

0 < y < 6, 



in which P (cos a) and Q (cos a) are the unknown spectral functions to be determined 
from the boundary conditions at y = 0 . These are 
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{Bl) Continuity of the tangential electric fields at y = 6. —oc < x < oc 
(B2) Continuity of the tangential magnetic fields at y = 8, x < 0 
(B3> Vanishing tangential electric field on the perfectly conducting half- 
plane at y = x > 0 


and we note that the boundary condition at the slab surface is implicitly taken into 
account by the representation (3.24) and (3.25). 

The application of (B1)-(B3) in conjunction with (3.24) and (3.25) is straight- 
forward. It results in the set of equations 


P (A) = —i'iQ (A) [l -v l R{\)e- j2k ° sV ^ 
r Q (A) v 2 (A) e~ jk ° xX d\ = 0 ; x < 0, 

J — oc 

r Q (A) t ’3 (A) [l - viR (A) e- 2jW ' /I=F ] e~ jk ° xX d\ = 

-f 3 (A 3 )/Mle^ [l - i h R( A„) e- w v^ 


(3.26) 

(3.27) 

gjfcox.Xo ;x > 0(3.28) 


where we have set A = cos a, X Q = cos <p Q and 


V2 (A) = 


1^3 (A) — 


1, E, polarization, 

1 / >/ 1 — A 2 , H z polarization. 
l/y/T— A 2 , E z polarization, 


(3.20) 

(3.30) 


l 1, H z polarization. 

These are sufficient to obtain a solution for C(A) and P( A). However, before pro- 
ceeding, it is necessary to rewrite certain terms in the integrands of (3.27) and (3.28) 
as products of '‘upper” and “lower” functions, that is, functions free of poles, zeros, 
and branch cuts in the upper and lower half A planes, respectively. In the process of 


doing so, we introduce the definitions 

O ■ 


'-’2 ( A ) d = v 2 (A) v 2 (A) 


(3.31) 
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in which 


MA) 


1 - i-iR (A) 


3 — 2 j k 0 6\ 


= t’ 3 + (A)^(A) 

dej Ltv (A; 6) (A; 6) 

L,(X)U,(X) 


(3.32) 

(3.33) 


£«, (A; 6) U w (A; 5) = .4 (A) (l + tq e~ 2 ^ s ^) 

+\/l-A 2 5 (A) f 1 - n e -2jWvCTj 
(A) n (A) = .4 (A) + Vl -X 2 B (A) . 


In these, L s (A), £u,.(A;<5), v 2 (A), v 3 (A) are lower functions while U s (A) ,U W (X: 8). 
vT (A), (A) denote upper functions. We also note that L W U W is a function charac- 

teiistic to the loaded parallel plate waveguide and its zeros correspond to modes in 
the \\a\eguide. On the other hand, L S 0 S is a function characteristic to the grounded 
slab with its zeros corresponding to the surface wave modes supported by the slab. 

Substituting (3.31 )-(3.35) into (3.27) and (3.28), we have 


/_>> 
/ 


Vo 


A) yj" (A) e~ jkoxX d\ = 0 ; x < 0, 


3.36) 


QW 4 W 4 (A) = 

L s ( A J U 5 (A) 


-I'ti-V 


(Ao) v /I~T eJ ».^ £1 Mhi 01 . x > 0 (3.37) 


These coupled dual integral equations can be solved for the unknown spectra by 
examining the analytic properties of the integrands, and the reader is referred to [10] 
and [58] for a more explicit description of this process. From (3.36) and (3.37) the 
unknown spectra are determined to be 


Q( A) 
P (A) 


MA 0 ;£)C,(A) 

2~jO w (A; <5) L 3 (A„) 


-i’iT a ,(A;^)L w (A 0 ;^) 

27tjL 3 (X)L,(X 0 ) 


4(A 0 )u 3 -(A 0 ) ^1-A 2 
v i (A) Vj (A) (A + A 0 ) 

4 l\U£(Kh/^ 

v 2 (A) Ug" (A) (A + A c ) 


Jk 0 8y/\-\l 


c jkoSWl-Xl 


( 3 . 38 ) 


(3.39) 



00 


Recognizing that Pdd{ A) = P{ A), we may then set <5 = 0 and substitute for the 
polarization dependent functions to obtain 


2njL a (A) L 3 (A 0 ) (A -t- A 0 ) 

L w (A) L w (A 0 ) y/l + Ay/1 -f A 0 
2rrjL s (X) L a (X 0 ) (A + A 0 ) 


(3.40) 

(3.41) 


where the superscripts e and h refer to the spectra associated with E, and H z po- 
larizations, respectively. 

To solve for the field coupled into the regio.. 0 < y < 8, x > 0, (3.21), (3.33) and 
(3.3S) are substituted into (3.24) to obtain the integral expression 

.+ i \ \ M 1 /l _ \ 2 T / \ c\ o i*o fi [\/ 1-A i-'/l- A2 j 


f: 


J ~> o 2") tC (A) v? 


vt(\ 0 )v;(\ 0 )Jl-\ 2 0 L w ( Xo-J) e 


jv.7 (A) vi (A) (A + X 0 ) L s (A„) I, (A) U w (A; 8) 

sin (kyyj 1 — A 2 ) , 

2j.4 (A) — = = = — - + 2 B (A) cos (k 0 y\fl - A 2 ) 


• jkx\ 


d\(3A2) 


for 0 < y < S. This can be evaluated by closing the path of integration via a semi- 
infinite contour in the lower half A plane. The sum of the residues of the captured 
poles then yields 


■ v f sin (k 0 yJ\ — A;:) / / 

F! = T C\ (A.) { JA(K) ; y + B (A„) cos ( k^ 1 - AJ 

1 y v 


n= 1 




(3.43) 


where {A n } are the zeros of U w (\]6), 


^ n “ c r / \ \ r / \ \ rr/ < \ . c\ „ - / \ \ _+ ( \ \ / \ t \ 


and 


L, (A n ) L, (A 0 ) U' w (A n ; 6) v 2 - (A n ) u+ (A n ) (A n + A 0 ) ’ 

(3.44) 


bl(X n ;S) = dUw 


d A 


A=A„ 


Substituting (3.29) and (3.30) into (3.44) with 8 = 0, we obtain the more explicit 
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forms for C,(AJ as 


cuk) = -- — ~h d±l 

L * L ° ( A °) twi'K) \ a + A n ' 

C n(^o) = IhdAA v/T+Xv/TTX 

lak) l s (A o) w w (a„) rrx; — 

for E s and ff. polarizations, respectively. 

3.3 Reflection and Launching of a Waveguide Mod. 

Consider "mv then n‘ A waveguide mode field (for 0 < y < 6) 


(3.47) 


F zn = 


sin / k 0 yy/l-\2 ) 

jA (An) L -7— +B{\ n ) cos (k 0 yyJT^T?) 


pj A n 


Fxn = [ l 'l- 4 (An) cos (kolfJl - A=) + j Vl B(\ \f \2 ■ (, r \1 

\ n) TJt-lO [Anhjl - A^Sin \k 0 yyj\ - A 2) C i*o*A n 

. .. L (3.48) 

incident at the waveguide opening (Figure 3 9 ( e n .. , . , 

° l i a ure «J-(e)). lhe radiated fields due to this 

excitation may he again represented by (3.24) - (3.25). Subsequent application of 

U ‘ e bOUn<iary C ° nditi0nS (B1WB3 > th ™ ‘b« -1-1 inte graI e, ualions (u , th the 

usual transformation to the A plane) 

Lw (A; S) U w (A; J) 

^ -e~ }k ° xX d\ = 0 ; z>0 

pjk 0 x A, 


VV 


here 


/ Q ( A ) Vj ( A ) v - ( A ) - w (A: £) V* 

Ls (A) U a (A) 

J_ x Q (A) f 2 + (A) y- (A) e-Fo^dX = - 

2 cos (M^/F- Ajj) 


; x < o. 


(3.49) 

(3.50) 


f A (A n ) , E z polarization 
r 4 = 4 


(3.51) 


, . ^£(A„) , H z polarization. v ™' 

he solution of those proceeds in a manner parallel to th* 

r parallel to the previous case. The result- 
ing spectra are determined to be 


0(A) = 


4t 


4 cos (Uv/T- Xl)L a (X n 


-v 4 LJ\ n ;S)UJX) v - (A n ) 


P (A) = 


)CC W (A; £) c 2 (A) y- (A n ) t-3 + (A) (A + A„] 
Uit’ 4 Z u; ( A n ; S)L W (A; A) (A„) 

(A) V 2 (K)vf (A) (A + An) 


cos (Wi - A5)i.(A„)i, (A) 


(3.53) 



with L w . L s , v ± . etc. as defined in (3.34), (3.35) and (3.31), (3.32). 


Substituting (3.4). (3.29). (3.30), (3.51) and (3.33) into (3.53) and setting <5 = 0 


we have 


Pn{ A) = 


-A(\ r 


L W (\)L W ( A n ) 


ef (A) = B(K )— f 


— A£ j X, (A) L a (X n ) A + A n 
L w { X) L w (\ n ) \/l 4- A\/l + A n 


(3.54) 


(3.55) 


^jL a (A) L a ( A n ) A + A n 
corresponding to the spectra for the E z and H z polarizations, respectively. The 
modal field reflected back into the guide may be computed by substituting (3.52). 
(3.21) and (3.33) into (3.24) and employing the usual transformation to the A plane 
to obtain (for 0 < y < 6) 

t ’4 U 3 (A n ) 


F! 


- L 


4 irj cos (M\/l ~ A*)u 2 (A) v 2 (A n ) t ’3 (A) (A + A n ) 
L w ( \ n ;6) 

L a (A) L a (A n ) U w {\\ 6) 

sin {k 0 yyfl - A 2 ) 


2jA(\) 


4 - '2B (A) cos (k 0 y\J 1 — A 2 ) 


- ]k ° xX d>( 3 . 56 ) 


Vl- A 2 

As in the case of coupling, this integral can again be evaluated by closing the path 
of integration in (3.56) via a semi-infinite contour in the lower half A plane to obtain 

sin 


FU = £ Rr 


m = l 


jA (A m ) 


fin (kpyy/l - A^) 


where R mn are the mode reflection coefficients given by 


+ B (A m ) cos (k 0 yyjl - A 2, 


> e 


“ J ko^\m 


(3.57 


1^3 (A n ) 


^ _ 

cos [k 0 8yjl - A^juT (A m )uT (A„)U 3 (A m )(A m + A n ) 


r w (An! ^) 


L a (A m ) L a (A„) U’ w (A m ; <5) 
When 6 is set to zero, this reduces to 

A (A n ) L w { A„) 


(3.58) 


/? 


e 

mn 


1 — A 2 L a (\ m )L a (\ n )U^(\ m ) A m + A„ 


3.50 
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Rln = B(\ n ) 


L w (A n ) 


\/l 4- A m \/1 + A n 


L s [\ m )L,(\ n )U' w {\ m ) 
lor the E, and II z polarizations, respectively. 


Am “1” A r 


(3.60) 


3.4 Computation of Spectra for Material Insert in a Per- 
fectly Conducting Ground Plane 


We now have all the necessary components required for constructing the spectra 
P's associated with e multilayer slab recessed in a ground plane as defined in 
(3.1 i' Substituting (3.40), (3.46), (3.54), (3.59) and (3.11) into (3.14) we obtain the 
E. polarization result 


P'A A- A.) 


where 


f,r(A)U,U yT^Ty/T^X 
L a ( A) L a (A 0 ) A + A 0 


1 + 


,v ,v 

EEr 

m = 1 7i=l A o 


V e 

mn 



(3.61) 


f' e _ ^ 777 1 Lw_ (A m) E ( E T 

mn yrrT 7 i,(A m )i s (A n) [/;(A n ) vl mV1 AnV 

Kn = {(/] + [Rmn]}^. 


(3.62) 


(3.63) 


For H z polarization, the spectra may be obtained by substituting (3.41), (3.47), 
(3.55), (3.60) and (3.11) into (3.14) to find 


PrA A,A 0 ) 


(A) Z , u , (A 0 ) \/l ~h A\/l + A 0 
-njL, (A) L, (A 0 ) A + A 0 


N N i ~rh 

m = l n=l p 



(3.64) 


where 


\/h R { \ \ ui ( A m ) 

mn " 1 m L a (A m ) L a (A n ) U' w (A„) 


\Jl + ~Kn\f 


i + KVL 


Vi = {[/] - [/?„„]} 


-1 

mn ' 


(3.65) 

(3.66) 


Expressions (3.61) and (3.64) can now be substituted into (3.13) and the resulting 
integral can be evaluated by the method of steepest descents to yield the far 
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non-uniform diffracted field 

F s (cos o. cos o a ) 


t~^ koP 

~ , / P TS (cos o. cos o, ' 


( 3 . 67 ) 


\ / » u * rs \ r '"U/ /— 

V VP 

where (p. o') denote the usual cylindrical coordinates. In (3.67), F s and P rs refer 
to Et and P* s in the case of E, incidence and to Z 0 Ht, P k s for H, polarization. 
Although not apparent, (3.67) is reciprocal with respect to cos d and cos do- as it 
should. We also note that P rs ( A, A r ,) is a combination of an inhomogeneous solution 
(direct diffracted term) and a sum of homogeneous solutions (modal contribution). 
It may also be easily shown from the asymptotic behavior of (3.61) and (3.64) that 
the homogeneous terms do not affect the edge condition. 


3.5 Specialization to the GIBC Representation 


To obtain numerical results, we must first provide expressions for the multilayered 
grounded slab reflection coefficient (i.e. -4(A) and B{ A)), as well as the associated 
split functions and corresponding complex roots (waveguide modes). To accomplish 
these tasks in a simple manner, we consider the general GIBC approximation to 
R( A). This amounts to setting 

.4(A) = ]T.4„(l-A ! )” 

n = 0 

Nq 

B (A) = £a,(l-A J )". (3.68) 

n = 0 

where A n and B n are constants specific to the multilayered slab and are given in 
chapter two. Introducing (3.68) into (3.21) yields 

- AT - 7I^En V fo^n(l ~ A 2 )" ' 

.En=0 A n (1 - AT + vT^FEnlo B n (1 - A*)". 


R( A) = — 


(3.69) 
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and note that for a given order of approximation. N A and Ng are finite and in general 
A a = As or A a = Afl + 1, with the order of the condition equal to 


A's = max {2N A ,2N g + T 


:3.70) 


To evaluate the split functions L W (X]6) and U w (A; 8), (3.68) is substituted into 
(3.34) and by setting 8 = 0 we obtain 

’~\A — A 2 ( 1 A 2 ) 1 ^polarization 

2 Ai U — A 2 ) n H z polarization. 

The split functions are then trivially obtained as 


L W {X)U W (\) = 


(3.71) 


t r w (\) = L w (-\) = 


f - B Sb \J 1 — A Iln=i \/l — sn — a] E z polarizati 




tion 


H z polarization, 

(3.72) 


where / m n/ 1 — < 0 and ((„} are the A A zeros of Hn=o B n \ n for E, polarization 

°r Hi i=o 3 n \ n tor H z polarization. It is apparent from (3.44) and (3.72) that the 


pertinent waveguide mode propagation constants are given by fc 0 \/l — Tr 
In a similar fashion, we may substitute (3.68) into (3.35) to obtain 


N, 


N B 


n—0 


n=0 

= Y:SnVT^ n 


2n+l 


n=0 


- *°n(i + AEZT 


(3.73) 

(3.74) 


w 


here 


n is even 


Sn = 


An / 2 

B( n - 1)/2 n is odd, 


Ns 


w Si 

{7n} = {zeros of the polynomial y — ( — 

/=o So 


(3.75) 
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The factorization of (3.74) is again trivial upon making use of the well known splitting 
germaine to the impedance half plane problem [34]. Noting that 


1 + t] v 1 — A ■ 


def 


yi-A : 


K+ (A: 7]) K- (A; 77) 


(3.76) 


we have. 

(vT^a)- Vj '/5o 

nS,A'+(A;l/7„) 

where 


r,(A) = 


(3.77) 


AV (A: t)) = A'_ (—A: ?;) 


A'+( A’?/) Ae(t?j>0 

|jA/ (A - ^'1 - 1/e) A' + (A; t]) | Re{y) < 0. 


and 


\J2 sin a/2,/^ 
v ' V v \ 

AP- ( 3 7T / 2 - a 

-9)*A*/2- 

- a + <?)]* 

«« (x/2) 

1 -f V 2 cos I 


1 + \/2 cos | 

f^)j 


In the above Re(rj > 0). Im [yjl - iff) < 0, 8 = sin _1 (i7),0 < Re{6), and is 
the Maliuzhinets function [22], whose evaluation in algebraic form has been given 


in [50]. Whereas the zeros of U u . (A) represent the waveguide modes, the zeros of 
L s ( A) (which are the poles of A'_) correspond to the surface waves supported by 
the material layer. Although not required in this analysis, these are easily extracted 
from (3.78) and (3.79). 

The expressions (3.68) through (3.79) provide a complete description of a GIBC 
implementation and permit the simulation of any multilayered coating. We remark 
that a unique GIBC modeling of a given coating does not exist; in fact one may 
employ GIBCs of substantially different character to simulate the same configura- 
tion. This point is discussed in the following section and some numerical results are 
provided for illustration purposes. 
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3.6 Numerical Results 

In this section GIBC simulations of various material inserts are presented and 
compared with exact results available for the case of a single laver. Due to its 
greater interest, data is presented only for the H z polarization case. The GIBC 
employed here-in are given in chapter two. These are valid for arbitrary multilayer 
coatings and are therefore suited for this application. In particular, these multilayer 
GIBC are synthesized by combining the component-layer GIBC in an appropriate 
manner, pointing to the necessity of understanding single layer simulations in order 
to construct multilayer ones. 

Figure 3.3 shows the far zone pattern of a single layer insert (e = 2 — ;.0001,/i = 
1.2.7 = .2A) modeled by various “low contrast” GIBC (i.e. those GIBC which 
improve as the layer thickness or index of refraction decreases). To illustrate their 
relathe contributions to the far zone pattern, the direct diffraction and modal contri- 
butions have been isolated in Figures 3.3(a) and 3.3(b), respectively, with the overall 
result presented in Figure 3.3(c). We note that for this low contrast GIBC, an 8 th 
order simulation provides a reasonable approximation to the diffraction pattern. 

In Figure 3.4, both the thickness and the index of refraction have been increased 
in a low contrast simulation of a (single layer) material insert with e = 3.5 — j.0001, 
ft = 2.0, and r = .4A. In contrast to the previous figure we now observe that a 
-0 order simulation is required to obtain a converged result. This degradation 
with increasing index of refraction proves typical of low contrast simulations and 
illusttates the need for other types of GIBC whose performance improves in this 
range of material parameters. 

In Figure 3.5 the same material insert corresponding to the data of Figure 3.4 
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(b) 




Angle In degrees 


(c) 



Angle In degrees 


Figure 3.3: H z polarization backscatter echo width for a material insert with r = .‘2A, 
e = 2 — j. 0001, fx = 1.2 modeled by low contrast GIBCs (see Table 3.1 for 
an explanation of the legend entries), (a) direct diffraction component 
(b) modal component (c) composite. 
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Igure 3.4: /^ polarization backscatter echo width for a material insert with r = 4A 
e - 3.5— 0001, n = 2 modeled by low contrast GIBCs (see Table 3 1 for 
an explanation of the legend entries), (a) direct diffraction component 
(b) modal component (c) composite. 
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is simulated with high contrast GIBC (i.e., a GIBC which improves as the index of 
refraction increases or as the layer thickness decreases). In contrast to the data in 
Figure 3.4. we now observe that only a second order high contrast GIBC simulation is 
required to accurately evaluate the far zone scattering. This difference in performance 
between the low and high contrast GIBC stems from the type of approximation 
employed in their derivation and the reader is referred to chapter two for a more in 
depth discussion. Some insight on the type of simulation provided by the low and 
high contrast GIBC may be gained through an examination of the waveguide modes 
predicted by the different simulations. These are presented in Tables 3.1 and 3.2 for 
the single layer simulations corresponding to the data in Figure 3.4 and 3.5. 

The exact modes given in the table are generated by the equation 

A n = Jtrfir- (-^) 2 , n =0,1,2,.... (3.80) 

We observe that as the order of the low contrast simulation is increased, the data in 
Table 3.1 reveal that the waveguide modes are “picked up 1 ’ in a sequential manner 
corresponding to increasing n in (3.80). On the other hand (see Table 3.2), the 
high contrast GIBCs pick up the n = 2 exact mode immediately and then “branch 
off'’ to pick up the other modes. The discrepency in pattern convergence between 
Figures 3.4 and 3.5 clearly suggests that the n = 2 mode is the most significant in 
terms of diffraction (for this particular configuration). We explain this physically by 
noting that the n — 2 mode may be resolved into its constituent rays which strike 
the interface at a characteristic angle (say which is greater than the critical 

angle 9 C of the material insert. On the other hand, the n = 0 and n — 1 modes are 
associated with characteristic angles less than the critical angle. This implies that 
upon coupling into the slab, the n — 2 waveguide mode is partially transmitted into 
free space while the lower order waveguides modes remain bound. 


c/A in dB ^composite) o/k in dB ^direct diffraction component) 
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Figure 3.5: H z polarization backscatter echo width for a material insert with r = ,4A, 
e = 3.5— j. 0001, [i = 2 modeled by high contrast GIBCs (see Table 3.2 for 
an explanation of the legend entries), (a) direct diffraction component 
(b) modal component (c) composite. 






Exact Modes 

L. C. ( 1 , 2 , 1 ) 

L.C.(5,6,3) 

L.C.(9,10,5) 

L.C.(13,14,7) 

L.C.(17,18,9) 

L.C. (21,22,1 1) 

2.6458 - jO 

2.6458 - jQ 

2.6458 - >0 

2.6458 - >0 

2.6158 - >0 

2.6458 - >0 

2.6458 - yO 

2.3318 ->0 


2.3324 ->.1518 

2.3303 - >0 

2.3318 - jO 

2.3318 — >0 

2.3318 - jO 

0.8660 - jO.OOOl 


-2.33 >4 - j. 151 7 

1.7608 - j.0001 

1.1583 - j.0001 

0.8782 ->.0001 

0.86612 — >.0001 

0 - >2.6575 



2.3156 — >.1518 

1.5228 - >1.4058 

0.8198 - >1.9935 

0.3937 - >2.6289 

0 — j 4. 24 64 



-2.3156 — 7-1517 

-1.5227 - >1.4058 

-0.8197 ->1.9934 

-0.3936 - >2.6289 

0 - >5.6624 




3.0782 - J2.1345 

2.2444 — >2.6318 

1.6576 -j 3. 1906 

0-;7.017S 




-3.0782 -> 2. 1344 

-2.2444 - j2.6318 

— 1.6576 — 73. 1906 

0 - jS.3404 





4.1090 - j2.8858 

3.1642 - j3.5300 

0 

1 

. 

w 

--j 





-4.1090 - >2.8858 

-3.1642 — >3.5300 

0 — 7 10.9348 






5.2491 -J3.5008 

0 - >12.2168 






-5.2491 - >3.5008 


Table 3.1: Low contrast approximation to waveguide inodes for a layer with e = 
3.5 - j.0001,/x = 2.0, and r = ,4A. For each low contrast boundary 
condition, the three numbers of the column headings indicate the ordei 
of the approximation in thickness r, the order of the resulting boundary 
condition, and the total number of modes (see chapter two). 



Exact Modes 

H. C.(0,2,l) 

H.C.(1,4,2) 

H.C.(2,6,3) 

H.C.(3,8,4) 

H.C.(4,10,5) 

2.6458 - >0 




7.9636 - >.0032 

3.2862 - >.0001 

2.3318 ->0 


5.5825 ->.0023 

2.2571 ->.0001 

2.2813 - >0 

2.2714 — >0 

0.8660 - >0.0001 

0.8667 - >.0001 

0.8649 - >.0001 

0.8660 - >.0001 

0.8660 - >.0001 

0.8660 — >.0001 

0 - >2.6575 



0- >2.1309 

^ 0 - >2.0694 

0.6002 - >2.2920 

0 - >4.2464 





-0.6002 - >2.2920 


Table 3.2: High contrast approximation to waveguide modes for a layer with e = 3.5— 
>.0001, /x = 2.0, and r = .4A. For each high contrast boundary condition, 
the three numbers of the column headings indicate the approximation in 
the index of refraction k - 1 , the order of the resulting boundary condition, 
and the total number of modes. 




69 


The above hypothesis may be tested by computing exact solutions in which the 
selection of the included waveguide modes parallels the order in which they are 
picked up depending on whether a low or high contrast GIBC simulation is employed. 
Figure 3.G depicts a high contrast simulation of a single layer having e = 11 — 
j.0001,^ = 7, and r = .4A. The 10-mode result is a pattern obtained by adding in 
modes sequentially as determined from (3.80), thus paralleling a low contrast mode 
selection scheme. On the other hand, the single mode result contains the contribution 
of only the n — 7 mode (the mode with 9 m > 0 C ), thus, paralleling the high contrast 
mode selection criteria. This clearly verifies that the most significant waveguide 
modes are those that are “visible’’, i.e. those with 9 m > 9 C . 

Finally, Figure 3.7 provides a simulation of a three layer insert composed of two *. 
mgh contrast layers (e = 11 — j.0001, = 7, and r = .4 A and e = 3.5— j. 0001, ji = 2.0, 

and r = .4A) placed beneath a low contrast layer with e = 2 — j.OOOl,^ = 1.2 and 
r = .2A. These are precisely the layers considered earlier in isolation. One might, 
therefore, expect that the order of a GIBC which provided a converged result for 
the single layer simulation will also provide an equally acceptable simulation when 
the slab is part of a multilayer stack. For the case at hand this is indeed true, as 
evidenced by the converged 9 th order result. We also remark that the presence of 
the two high contrast layers enhances the modal contribution to the total diffraction 
when compared with the single layer data given in Figures 3.3(b) and 3.3(c). 

3.7 Summary 

In summary, the scattering from a vertically inhomogeneous slab recessed in a 
ground plane was obtained through application of the generalized scattering matrix 
technique in conjunction with the dual integral equation approach. The solution was 
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Figure 3.6: H z polarization backscatter echo width for a material insert with r = .4A. 

t = 11. - j.0001, fi = 7 modeled by high contrast GIBCs (see Table 3.2 
for an explanation of the legend entries). 
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Figure 3.7: H z polarization backscatter echo width for a three-layer material insert 
with (r! = .4A, ei = ll.— j.0001, fii = 7), (t 2 = .4A, e 2 = 3.5 — j.0001, 
H 2 = 2), (r 3 = .2A, t 3 = 2. — j.OOOl, /z 3 = 1.2). In the legend entry 
a.b,c(d,e), a and b denote the approximation in /c” 1 of the high constrast 
layers 1 and 2, c denotes the approximation in r of the low constast layer 
3. while d is the order of the composite GIBC and e is the total number 
of modes. 







72 


specialized to the case of a multilayered slab simulated with a generalized impedance 
boundary condition (GIBCj. Results were given for various single layer inserts and 
it was seen that in the case of materials having sufficiently high index of refraction, 
high contrast GIBC simulations converged more rapidly (with respect to the order 
ot the GIBC) and performed better than low contrast simulations. Finally, results 
were presented for a 1A thick lossless three-layer insert containing both high and low 
contrast layers. It was shown that the simulation converged at the point predicted 
by the individual layer simulation suggesting a method for constructing multilayer 


simulations. 



CHAPTER IV 


GENERAL SOLUTION OF THE 
DIFFRACTION BY A MATERIAL 
DISCONTINUITY IN A THICK 
DLELECTRIC/FF TRITE SLAB 


In this chapter and the following one it is demonstrated that a GIBC/GSTC 
sheet characterization can yield a unique solution when supplemented with certain 
conditions at the sheet discontinuity which do not require an apriori knowledge of 
the edge fields. As a vehicle for presenting this solution procedure we employ the 
dual integral equation method to consider the plane wave diffraction by a discon- 
tinuous distributed sheet (see Figure 4.1(b)). This very general model is capable of 
representing material half-planes, material junctions, and material discontinuities on 
grounded structures, such as those shown in Figure 4.2. In addition, a distributed 
sheet model typically renders the same degree of accuracy as the usual infinitely-thin 
sheet, but with a lower order condition. It is, therefore, of much practical interest. 

In the first section of the chapter, the GSTC representation of the distributed 
sheet discontinuity is used to develop dual integral equations in terms of the un- 
known spectral functions proportional to the sheet currents. These equations are 
then solved in the standard manner to yield expressions for the spectral functions in 
terms of unknown constants. The constants are dependent on the material and geo- 
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Figure 4.1: (a) Distributed sheet, (b) Distributed sheet discontinuity. 




Figure 4.2: Geometries modeled by a discontinuous distributed sheet, (a) Material 
half-plane, (b) Material-material join, (c) Grounded join. 
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metrical properties of the discontinuity and their specific value is identified here for 
several discontinuities by comparison with a few known diffraction solutions. This 
demonstrates the validity of the presented solution, but in general, the determination 
ot the constants requires the enforcement of additional constraints demanding field 
continuity accross the layer discontinuity. The development of these conditions and 
their use in solving for the constants are presented in chapter five. 

4.1 Dual Integral Equation Formulation 


Consider a distributed sheet of thickness r illuminated by the plane wave 


Fine — ^ 


jk(z cos o 0 + y sint^o) 


^2,tnCl E z polarization, 
Z c,H z i nc , H z polarization, 


(4.1) 


as shown in Figure 4.1(a). The excitation (4.1) induces reflected and transmitted 
fields which are explicitly given by the properties of the distributed sheet. If this 
sheet models a symmetric slab, then an appropriate GSTC representation is formally 
given by (see chapter 2) 


) { F+ ~ F } + (~Tr) { dy [ F+ + F ] } = 0, -oo < X < oc 

(4.2) 

(-ir) + F '} + fa (-ir) [f + -f-}} = o, 

(4-3) 

in which F is the total field, F ± = F (x, y = ±r/2), dxF ± = f^F (x, y = ±r/2), and 
dyF- = j^F(x,y) | y= ± T / 2 . Also, U} 2 (-^) are differential operators which operate 
on the field quantity in the curly brackets, and are finite polynomials in whose 

coefficients depend on the slab modeled by the distributed sheet. To maintain the 
generality of the solution, the U}- operators are left in symbolic form and the reader 
is referred to chapter two for their explicit representation in terms of the material 



constants and thickness of the layers comprising the modeled slab. In general, the 
order of U xx (i.e. the highest derivative present) is usually the same or one more than 
that of l( x2 and similarly the order of U 2 1 is the same or one more than the order of 
l.{ 22 . Thus, we may define the orders of the GSTCs in (4.3) to be 

X° dd = maximum j order of U{ x (A 2 ) . 1 + order of U\ 2 (A 2 ) } 
yeten _ maximum | order of Z /21 (^ 2 ) ■ 1 + order of U\ 2 (A*)| 

(4.4) 


The reflected and transmitted fields may now be easily determined by employing 
(4.3) to find 

F refl = fl ie JM*co S 0 o -ysin* o ) (4.5) 

p __ p gjk[x cos + y sin 0o) (4.6) 


in which Ri and T\ are the reflection and transmission coefficients, respectively, and 
are given as 


pjkr sin 4> 0 

R x = — - — [ R\ ven + R° x dd 

pjkr sin <b 0 . 

p y pcvcn R 0( *d 


(4.7) 

(4.8) 


with 


sin <pJJ 22 (cos 2 <f> 0 ) — U\ x (cos 2 d> 0 ) 
sin <j> JA 22 (cos 2 4> 0 ) + W 21 (cos 2 <p Q ) 
sin <t>Mn (cos 2 <f> 0 ) - U\ x (cos 2 <ft Q ) 
sin 4>Mi2 (cos 2 <p Q ) + U\ x (cos 2 4> 0 ) 


(4.9) 

(4.10) 


We remark that in (4.9) and (4.10), U} } (cos 2 d> 0 ) now represent simple polynomial 
functions in cos 2 O a . since —dx 2 /k 2 = cos 2 <p Q in view of the field expressions f 4.5 1 


Consider now the case where the right half of the distributed sheet in Figure 4.1(a) 
is replaced by another sheet of the same thickness, but of different properties, as 
illustrated in Figure 4.1(b). The GSTC representation of this modified sheet is 



for — oc < x < 0 and 


{f + + F"} + {dy [f + - F“] } 


0 

0, (4.11) 


W?1 

u\ x 





0 

o, (4.12) 


tor 0 < x < oc, where the superscripts 1 and 2 distinguish the left- and right-hand 
sheets, respectively. Referring to our previous discussion, the orders of the right hand 
side GSTCs are given as 


*\ 2 dd — max { order of U\ x (A 2 ) in A, 1 + order of U\ 2 (A 2 ) in A j 
* v r en = max { order of U\ x (A 2 ) in A, 1 + order of U\ 2 (A 2 ) in A} . (4.13) 


The modified right hand side sheet induces a scattered field F, in the presence of 
the excitation (4.1), and the total field can now be represented as 

Fmc + Frefl + F, y > t/2 

Ftran + F, y < r/2 

where F s is the unknown scattered field in the region \y \ > t /2 and can be expressed 


F = 


(4.14) 


as either [10] [5 


F : 


i*,y) = f c 


— Podd (cos a) + P even (cos a) 


e -jfcsma(|j/|-T/2) e -;*;xcosaf^ a / j j-j 
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(b) 


Illustration of C contour in the complex a plane and (b) complex A 

e. 
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or 


F s ix.ij) = 


jPoddW + P tven ( A) 


- |y | -r/2) e - j kx A 


(4.16) 


V 1 — A 2 

upon invoking the transformation A = cos a. The contour C in (4.15) is the inte- 


giation path in the complex a plane shown in Figure 4.3(a), and its counterpart for 
the A plane is illustrated in Figure 4.3(b). Also, the spectral functions P odd ( A) and 
P‘ien ( A ) are directly related to the Fourier transforms of the unknown equivalent 
currents 


J odd — P (4.17) 

J even = F+ + P s . (4. IS) 

via the relations 

/ oo j \ 

Podd{ \)e-P* x ^== (4.19) 

vl — A* 

/ OO r ] \ 

Peven( (4.20) 
-oo vl — A 2 

Assuming that 

Jodi{x) ~ x 3 ° dd as x — > 0 

Jeven{x) ~ X 3even as X — ► 0 (4.21) 

with 0 < s odd < 1 and 0 < s even < 1, from (4.19) - (4.20) and the Abelian theorem 
we have 

Podd{ A) ~ \~ 3 ° dd as j A | — ► oo 
Peven (A) ~ A _, '’'«" as |A| -» oo (4.22) 

From the asymptotic behavior of P 0 d d { A) and P even (A), it may be easily seen that 
the integral representation (4.16) is well-behaved and convergent for all x and y. 
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Substituting (4.1), (4.5), (4.6), (4.14) and (4.16) into the transition conditions (4.11) 
and (4.12) we obtain 

— -fi-p = 0 (4,23) 

/” ff ""( v ) p '~” ( A ) e ‘ , ‘ lJ = °’ ( 4 ' 24 > 

for x < 0 and 


/“« 251“" (A 2 ) />,«„ (A) 


for x > 0. where A 0 = cos <? 0 and 


2sin4,e2* lX °^* r / 2lin *°Z 0<(t ((A|) 

S^A 2 ) 

2 sin 0 o e> klX ° e J kT / 2 * in ( A‘ ) 

^r* n (A^ 


(4.25) 

(4.26) 


Q\ dd (A 2 ) = U l u (A 2 ) + a/1 -A hl\ 2 (A 2 ) (4.27) 

geven = ^ ( A 2 ) + ^l-A^/ (A 2 ) (4.28) 

£f J (A 2 ) = w 2 (A 2 ) + >/T^W 2 2 (a 2 ) (4.29) 

graven ^ = ^ ( A 2 ) + (A 2 ) (4.30) 

Z 0<W (A 2 ) = [^li (a 2 ) ^12 (a 2 ) — ^12 (Aq) (Aq)] (4.31) 

Z e „ en (A 2 ) = [u x 2l (Ao) ^22 M - U\ 2 (A 2 ) ^ (A 2 )] . (4.32) 


Xote that strictly speaking, the integrals in (2.26)-(2.31) are not convergent because 
of the polynomial order of the integrands. This difficulty is common in analytical 
GIBC solutions and may be remedied by working with integrated field quantities 
as discussed in [39] and [44]. It has been shown, though, that the final solution is 
the same regardless of this remedy and for the sake of simplicity we will proceed 
with the solution of the dual integral equations as if all Fourier inverse and forward 
transforms existed in the classical sense. 

Equations (4.23) with (4.25) and (4.24) with (4.26) form two uncoupled sets of 
integral equations which are sufficient to yield a solution for the unknown spectra 


Poddi an d P even (A). Clearly, because of the similarity between the two sets of 
equations, once a solution for P^dd{X) is found, the corresponding one for P even (A) 
follows by inspection. 

4.2 Solution of the Dual Integral Equations 


In proceeding with the solution of P odd (A) we first rewrite (4.25) in a more suitable 
manner, viz. 



2<rj'Cf“(A„)(A + A.) / 


= 0 


(4.33) 


for x > 0. Closing now the integration paths in (4.23) and (4.33) by semi-infinite 
contours in the upper- and lower-half A planes, respectively, leads to the deduction 
that 


Gt dd ('\ 2 )P«u( A) 
vIPP 


U 0 ( A) (4.34) 

E 0< td (A) . ... 

2*j6i“W)(.X + K) A„) ° ’ 


where U 0 ( A ) and L 0 (A) are unknown functions regular in the upper- or lower-half 
of the A plane. Also, E odd (A) is an unknown entire function to be determined along 
with P od d{ A), f/ 0 (A) and L 0 ( A). To solve for these, it is necessary to exploit the 
analyticity properties of (4.34) and (4.35) in the different regions of the complex A 
plane. An important part of this process is the factorization of the Q functions in 
(4.34) and (4.35) into a product of upper and lower functions (that is, functions free 
of poles, zeros, and branch cuts in the upper- or lower-half A plane, respectively). 
This task is described in the appendix. We have 


Qi dd ( A 2 ) = G° l d + d (X)g o 1 d _ d (\) 

Gl vtn (a 2 ) = (A) Q[ v J n (A) 
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G° 2 dd (A 2 ) = 0 2 o "(A)£|"(A) 

Gr n (A 2 ) = Gll' n (\)G' 2 l' n (\) (4.36) 


where the snperscipts ± denote upper or lower functions, respectively. Combining 
(4.34) with (4.35) in conjunction with (4.36) and rearranging terms we obtain 

gf?(A)(7o(A) 5 (X)G it d ( A ) _ MA)g°l'(A) 4:3 - , 

Gfl d { A) (A + A.) G° 2 t d ( A) GfJ (A) 

where 

sm0 0 e^ 2sm ^Z odd (^ 0 ) E odd (X) ,, 

- (A) = wFw) E M (-\ 0 y (4 ' 3b) 

The second term of (4.37) may be easily split into a sum of upper and lower functions, 
and when this is done, (4.37) can be rewritten as 


gff (A) V 0 (A) 5(A)gff (A 0 ) = Lp (A) G±^_ (A) 

Gti d (A) + (A + A 0 ) Gt d (Ao) Gt- (A) 

S(X) r gfgjA) _ gff (A„) 1 
(A + A 0 ) [G 2 i d (A) (?2+ d (A 0 ) j ” 

The left hand side of (4.39) is now r regular in the upper half of the A plane while its 
right hand side is regular in the lower half of the A plane. By Liouville’s theorem, both 
sides must then be equal to a polynomial, and to determine the order of this poly- 
nomial it is necessary to examine the asymptotic behavior of the individual terms in 
(4.39). From (4.4) and (4.13) G 2 ± (A) ~ \ N ? dd / 2 , G{± (A) ~ \ N ° dd /2 , implying that the 
left hand side of (4.39) behaves as |A |( A C dd+A 2 <ili )/ 2-1 “ J °<M when |A| — ► oo, provided that 
E odd (A) behaves no worse than | A|6XLd +A odd)/ 2-, <> d <* as |A| — ► oo. The right hand side 
of (4.39) will then behave as the greater of |A| Ar > od<f-1-a ° <<d and | y | ( + N° dd )/2-i-s <Kid 

when j A | — * oo. In accordance with Liouville’s theorem, both sides of (4.39) will 

f N odd 4- N odd 1 T . . 

then be equal to a polynomial of order int < — 1 — ^ 1 — In terms of this 

unknown polynomial, we may solve first for U 0 (A) and subsequently for P odd ( A) to 
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find that 


Po4d( A) = 


J sindov/l-A" 


c jkrf2 sin <$> 0 


- 7r A + a, gfj (A) gtt d (\ 0 ) Qii d i\) Qi d + d {\ 0 ) 

/ \ jj { A ^ odd 1 “7 TI 

Z °“( X D r _A + E E «™.(A + A„)"(AA„r (4.40) 

■ c '0(ii i *o) m = 1 n=0 


In this, .\ 0 dd = int {1/2 (A*^ + A ^ dd 4- 1)}, and a mn are arbitrary constants as yet 
undetermined, and correspond to the coefficients of the polynomial resulting from 
the application of Liouville’s theorem. The chosen symmetric form of this polynomial 
is not unique but will be found most useful later in constructing a reciprocal form 
for P odd (A). 

Following a similar procedure we also obtain P even ( A) as 
_ . . 1 sin d>„\/l — A 2 p jkT/2sm<t> 0 


p , , . _ j sindovl— A 2 gjfcr/ 2sin0 o 

27 a + a 0 (A) g?™ (A 0 ) gz™ (A) ggr (J7) 

, . (A) v«n ~ 1 J^etun — 1 — m j 

Z e ven (^ 0 j — + ^ ^ b mn (\ + \ 0 ) m (\ \ 0 ) U ( 4 . 41 ) 

^evenK *o) m = 1 n=0 

with E even (A), A' et . en and b mn being the counterparts of E odd { A), N odd and a mn , 
respectively. We note that (4.40) and (4.41) imply that the powers s odd and s et , en 
governing the behavior of the equivalent currents J odd ( x) and J even (x), respectively, 


Sodd — 


1 N odd + N odd is even 

1/2 N odd + N odd is odd 

1 N even + N Ln is eve n 


(4.42) 

(4.43) 


( !/ 2 N Len + is odd 

To determine the unknown entire functions E odd ( A) and E even (A), we observe 
that the spectra P odd (A) and P even (A) must exhibit a reciprocal form, which may be 
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achieved by setting E 0 dd{ A) and E eve n (A) equal to any of the following functions: 

( Z odd ( AA 0 ) or 


E<m( A) = 


Z odd { A ) or 


(4.44) 


l^(A) 

r* even ( — AA 0 ) or 


E eve n(A) = < 


(4.45) 


Z even ( A ) ° r 

^ Zeven (A) 

where Z+ ddeven (\) and Z~ ddeven (A) are upper and lower functions satisfying the re- 
lation 

Z odd, even (A*) - Z^ ddeven (A) Z~ dd even (A) (-1.46) 

Taking into account the choices (4.44) and (4.45), we may substitute (4.40) and 
(4.41) into (4.15) and subsequently perform a steepest descent path evaluation to 
obtain for p — * oo (all surface wave contributions are neglected in this evaluation) 

e-jkp 


F (p, <p) ~ [D od d (d>, 4>o) + D even (p, <£<>)] 


y/kp/ 2tt 


(4.41 


where (/?, 0) are the usual cylindrical coordinates and jDodd(0, 0 O ) + £>even 0o) is 
the far zone diffraction coefficient symmetric with respect to 0 and <fi 0 . We have 

e -W4 sin0 o sin0 


D 0 dd ( 0 ? <Po) — 


27T cos 0 + COS 
e j/pr/2(sin^ 0 +|»in^|) 


G{t d (cos 0) £{1* (cos <p 0 ) (cos <p) Gt¥ (cos 4> 0 ) 

^ N ydd. — 1 N od d - 1 -m 

2^ (cos d>, COS </> 0 ) + ^ E Qmn ( COS ^ + cos ( COS ^ cos 0°)" 

m=l n=0 


(4.48) 


Deven (0? 0o) — 


e Jir / 4 sin0 o |sin0| 
2n cos 0 + cos 0 O 

p jkr /2{sm <}>o+\ sin 4>|) 


^even ( CQS ^ Qeven ( CQS ^even ( cos ^ Qcven ( CQS 

^ ^«vcn~ 1 ^euen 1 — fU 

Z even (cos <2>,COS d> 0 ) + Y1 b mn (cos<j>+cos<j> 0 ) m (cos<f>cos<p 0 ) n 

m = l n=0 


(4.49) 
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in which the functions Z 0 dd,even (cos o. cos <p 0 ) are given by (see (4.44) and (4.45)) 

Zodd (— cos <t> cos 4> 0 ) or 

2 odd (cos <p) Z~ dd (cos <p 0 ) or (4.50) 

2 odd ( c °s <f>) Z+ d (cos 0 O ) 

Z even (— cos <?cos <p 0 ) or 

Cn ( COS $) Z even ( COS 4>o) Ot (4.51) 

Z tven ( COS <t>) Z tven ( C0S <f>o) 

Because the above three choices for Z 0< id and Z evcn differ only by terms of the form 
(cos n - r cos <t> 0 ) m (cos <p cos 0 o ) n , it is immaterial which of them we choose, although 
one of the choices may likely lead to a more compact representation. Nevertheless, 
regardless of the choice of Z 0 dd and Z evcn , one is still faced with the determination of 
the unknown constants a mn a*nd b mn in (4.4S) and (4.49), respectively. These are a - 
manifestation of the non-uniqueness of the finite-order GSTC sheet model employed 
herein, and their explicit determination requires the introduction of additional con- 
straints pertaining to the physics of the problem as discussed in chapter five. In some 
cases, however, these constants can be determined by comparison with alternative 
diffraction solutions, and this is considered next. 

4.3 Discussion of the Solution and Some Applications 

4.3.1 Diffraction by thin single layer discontinuous slabs 

The diffraction coefficient given by (4.48) and (4.49) is very general and can 
model a wide variety of geometries. To check its validity, display its versatility, and 
assess the relative importance of the unknown constants, we consider several simple 
configui at ions which can be modelled with the proposed GSTCs. Their geometries 
are shown in Figures 4.4(a)-4.4(c) and include the single layer join, the material-metal 
join, and the material half plane, all of thickness 2 w. Herewith, these are modeled 
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Figure 4.4: Thin slab structures and their distributed sheet representation (a) 
Material-material junction, (b) Material-metallic join, (c) Material half- 
plane. (d) Equivalent distributed sheet. 



ss 


by a distributed sheet discontinuity of thickness 2 (w - w a ) (see Figure 4.4(d)) and 
although, in general, it is not necessary to employ a sheet thickness different from 
that of the slab, this is useful here for comparison with previous results. 

If the left hand side of the slab is assumed to be thin and having a low index of 
refraction, it may be modeled by a low contrast GSTC sheet. Thus, an 0(w 1 ,w]) 
approximation with terms of 0(w s w) neglected is sufficient for the representation of 
the operators or polynomials U\ y In particular, we have 

U\ x (- dx 2 jk 2 ) = 1 

^12 [—dx‘/k 2 ^ = jk(u\w — w a ) 

Mi (-*'V* S ) = Jk (^ - w.) + { (^ - w.) dx’ 

U] 2 ^—dx 2 /k‘'j = 1 (-1.52) 

where £i and fi\ are the relative permittivity and permeability of the left hand slab, 
respectively, and 

{ Hi, E z polarization 

(4.53) 

ci H z polarization 

Also, when w, = w, these are simply the transition conditions derived first by We- 
instein [59] and later by Senior and Volakis [39]. The corresponding polynomials to 
be employed in (4.27) - (4.32) are given by 


U\ x (— COS 0COS 4>o) 
U\ 2 ( — COS <f> COS 0 O ) 
U\ x (— cos 0 COS 0 O ) 

u 

>2 (— COS 0 COS 4>o) 


1 

jk (u x w - w a ) 

( weiHi \ ., fw \ 

jk y “ W a J + jk ^ W a J COS 0 COS 0 O 

1 (4.54) 


Incorporating these into (4.48) and (4.49) and setting 


Z 0 dd (cos 0, COS 0 O ) 


Zodd ( — COS COS 8>o ) 


(4.55) 



89 


Z even (cos 0, cos <p 0 ) = Z even ( - cos 0 cos <f> 0 ) 


(4.56) 


yields 


D odd (4>Z6 0 ) = - e . W4 Smg>oSm ° e jfcr/2(»«».+|,in»|) 

2jt cos0 + cos<p o 

7/ 2 2 ( — cos 0 cos 0 O ) — jk (u\W — tUj) Ll\ (— cos 0 cos 


D, 


U\ 2 ( — cos 6 cos 0 O ) — jk(u\W — w 3 )U xl (— cos 0 cos do) 
' M_ (cos 0; 7 0 ^- 1 ) Ml (cos 0 o ; 7 0 "’ 1 ) (cos ^ gTTd (cos 

S “ W4 Sitl 0ol Sin ^L.ifcT/2(sin0 o +| S in<(,|) 


(4.5' 


tven ( Qi Oo) 


2 7T COSO 
[a t + ao cos o cos 0 0 


S ^ n c ,?'fcr/2(sin tf> 0 + | sin <3 


k 03 [ilmsti d/_ (cos 0; 7m en ’ 1 ) A/- (cos 0, 


s 0 O ) 


+ 


cos Oo 

>] ^22 ( — COS 0 COS 0o) — U\ x (— COS 0 COS ■ 

*) A/_ (cos 0 O ; 7m en ’ 1 )] & 2 + ” ( COS 0) £ 2 + n ( cos 9 ° 

6jo (cos 0 4- cos 0 O ) ^ 

7 1 \ 


O' 3 n 2 m= lM- (cos 0;7m en,1 )iV/_ (cOS 0 O ; 7m' ' n,1 )j <?2+ en ( cos 0)£!+' " ( COS 0o) 


where the split function M_ (cos<£;7) is given in the appendix. 


<*x = 
a 2 = 

o=3 = 


(z ~ w ‘) 


jkw 

u 1 


( e iMi — 1) 


and 


(4.58) 


(4.59) 


odd, 1 


7 


even,! 

7l,2 


-J 


A: (u x w — in,) 

til ± \Jv\ + ikho (ei fix — 1) (to — Ujto,) 


(4.60) 


2 jk ( w — w 3 u x ) 

with 'y odd ° T even representing possible surface wave poles. To complete the definition 
of (4.57) and (4.58), the functions associated with the right hand side properties ot 
the slab (i.e. those functions with the superscript 2) must be specified. Referring 
to the configurations in Figure 4.4(a)-4.4(c), Table 4.1 and Table 4.2 provide ex- 
plicit expressions for the functions ZY 2 (— cos 0cos 0 O ), (cos 0) £f+ d (cos o 0 ) and 



Geometry, x > 0 

U\\ ( — cos <f> cos 4> 0 ) Ui 2 ( - cos $ cos <f> 0 ) 

&2+ d ( COS <f>) 02? (COS 4> 0 ) 

syOdd, 2 

Low Contrast 0(w,w t ) 

1 

jk(u 2 w - w t ) 

M + (cos <f>] 7 r 2 ) M+ (cos 4> 0 \ 7 , 0<W ' 2 ) 

-1 

Free Space Limit 

1 

jk(w-w,) 

M + (cos <t>] 7 1 oJd ’ 2 ) M + (cos <p 0 \ 7 ^' 2 ) 

Z 1 

k(w-w 9 ) 

PEC(Ez-pol) or PMC(Hz-pol) 

1 

0 

1 


PMC(Ez-pol) or PEC(Hz-pol) 

0 

1 

n/ 1 — cos — cos <j> 0 

■ - 


Table 4.1: Odd symmetry parameters for distributed sheet similation of right hand 
side material x > 0. See Appendix 1 for definition of A J + split functions. 




Geometry, x > 0 

(-cast cos fa), W| 2 (-cos^cos^ 0 ) 

a?r(cos^)air(«>^o) 

Low Contrast, 0(tu T u»,) 

M* = jk (“jf 1 - «’•) + 7 k cos 4> cos <j> 0 - w.) 

= 1 

** (c 2 /0 - l)nL, Ah (oos^; 7 r u>J ) Ah (cos^sr-*) 

euen,2 ± \J » Ak 7 w(i2»2 - 1 )(u/-*2 »"* ) 


Wltn 2 — — u;* uj ) 

Free Space Limit 

U\ x — jk(w - tA>.) ( 1 + COS <t> COS <t> a ) 

U\, = 1 

\/l - cos ^v/1 - cos (cos <t > ; 7 J' VCT ‘' 2 ) A/ + (cos <£ c ; 7 [ l ' cn * 2 ) 

with 'v tutn ' 2 — 

PEC(Ez-pol) (6 10 = 0) 
or PMC(Hz-pol) (6,o = 0) 

i/ 2 , = 1 

*4 = 0 

1 

PMC(Ez-pol) or 
PEC(Hz-pol) 

U\ x = 0 

*4 = 1 

\/l “ cos ^\/l "" COS to 


Table 4.2: Even symmetry parameters for distributed sheet similation of right hand 
side material x > 0. See Appendix 1 for definition of A/+ split functions. 
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tf 2+ (co !0)6 m (cos«„) terms. By edge condition considerations, all of the con- 
slants and have been set to zero except * 10 appear, ng in the definition of 

which is non-zero unless the right hand side slab is a PEC/PMC under an 
Hz/ H z excitation (see Table 4.2). 

4.3.2 Diffraction by a resistive-resistive junction 

When w = u„ and the material parameters of the slab geometries in Fig- 
ures 4.4(a) and 4.4(c) take the limits 


e l,2 


oo 


k 1,2 


jkw (e 12 - 1) 


oo 


0 


* 7 1,2 


jkw(n i,2-l) — ► — 

7 „m 


V?,2 


(4.61 ; 


the resulting configuration corresponds to coincident resistive and conductive sheet 
junctions [37], where ,f, 2 and denote the resistivities and conductivities of the 
respective sheets. In this case, the constant 4,„ is forced to zero by the edge condition, 
and (4.5/) and (4.58) simplify to 


D 0 dd (COS 0 , COS (j ) 0 ) 


sin d> d (cos d>, cos <^ 0 ) (l/r?? 

M ~ ( cos ft l hi) A/- (cos 0 O ; 1 /t?P) M + (cos 

D even (cOS <j), COS 0 q} = 


v^72 

(cos 0; l/T)?) M + (c 


tor £\. polarizat 



(cos <f>, cos d> 0 ) = 
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sin d d (cos 4 >, cos <t> 0 ) (l/r/| — I/77*) 

A/_ (cos 6 ; I/77*) A/_ (cos <? 0 ; 1 / ) A/+ (cos <£>; L/77I) A/_ (cos <? 0 ; 1 / 77I ) 


(4.64) 


£h Len (COS O. COS (p 0 ) - 


| sin <*>\ d (cos 6. cos o 0 ) {vT ~ VT) 

M- (cos <p\ l/r/™) M- (cos <p 0 ; l/j/j 71 ) A/ + (cos <p; l/r}™) A/+ (cos o 0 \ 


T7^ 3) 


for H. polarization, where 


d (cos <p, cos 4> 0 ) = — 


e • J ’ r/f4 sin <t> 0 
2 x cos <p + cos <p 0 


(4.66) 


Note that if both 77 *, t\\ -+ 00 (i.e., the resistive sheets disappear), then 
-Dodd (cos 0 , cos 0 o) and D ei)en (cos <?, cos <b 0 ) tend to zero for H z and D, polarizations, 
respectively. On the other hand, if both r/j 71 , rj™ — >- 00 (i.e., the conductive sheets 
disappear), then D 0£ fd (cos <j>, cos o 0 ) and D even (cos cos 4> 0 ) tend to zero for E z and 
H z polarizations, respectively. This is, of course, because resistive and conductive 
sheets scatter independently of each other. Consequently, the field diffracted by a 
resistive to conductive sheet junction is the superposition of the individual sheet 
contributions [39]. 


4.3.3 Diffraction by grounded metal-dielectric junctions 

Of the geometries shown in Figure 4.4, the diffracted field associated with the 
metal-dielectric junction (Figure 4.4(b)) is given in chapter three and can therefore 
be used to partially validate the derived solution. However, in order to study only 
the effect of the constant 6 10 , we need to exclude the odd-symmetry portion of the 
metal-dielectric join diffraction coefficient. To this end, we focus on the recessed slab 
geometry of Figure 4 . 5 (a), whose H z polarization diffraction coefficient is related 
(through image theory) to that of the metal-dielectric join by 


D t5 (cos <p, cos <t> 0 ) = 2 D even (cos 4>, cos <j> 0 ) 


(4.67) 
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with D ewm given by (4.58) along with the PEC(tf,-pol.) entries of Table 4.2. As 
noted above, for a right hand side PEC slab with E. polarization, the edge condition 
demands b 10 = 0 in this case. 

For the recessed slab geometry illuminated with an H z polarized plane wave, the 
GSTC (4.11) and (4.12) become the GIBC 



for — oc < x < 0 and 



for 0 < x < co, with the U\ : operators given in (4.52). Clearly, these GIBC can 
represent any of the configurations displayed in Figure 4.5(a)-(c) without regard to 
whether a stub (PEC or PMC) or not is placed at the junction. This information 
can only be carried by the constant 6io as the term distinguishing the diffraction 
coefficients among the geometries of Figure 4.5(a)-(c). Thus, the determination of 
bio must somehow involve the properties of the junction across its thickness and 
this is discussed in chapter five. However, since the solution of the configuration in 
Figure 4.5(a)-(c) are already available in chapter three, b \ o can be identified for each 
geometry by comparing (4.58) with the appropriate solutions given in chapter three. 
Upon setting w a = 0, we find 


ino stub 


'10 


• pec 
°io 


*r c 


= jkwJ— 


stub 


stub 


jkwJf 


j kw 
4c 


f (v^ - i) [m_ (v^r*”' 1 ) m - (v^r^r"' 1 )] + 1/2 


jkwJf 


jkw 

4c 


f (v^jm - 0 [m- m- f^Kor”' 1 )] - 1/2 

(4.70) 
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corresponding to the constants associated with the diffraction coefficients for the 
geometries in Figure 4.5(b), Figure 4.5(a), and Figure 4.5(c), respectively. We note 
that these constants are based on the choice of Z fven (cos 0,cos 0 O ) as given in (4.56). 
Had this function been chosen as 

Z even (cos 0, cos 0 O ) = Z~ en (cos 0) Z~ en (cos 0 O ) (4.71) 

the resulting constants would have been 

mo stub _ n 

°io — u 

,pec stub /Mi 

6 10 = -jkwJ — 

+ J kw \f^ 

ikf (v /e iA £ i - l) (M- ( N /ciMi)7i l ' eT1,1 ) M- (v^T,7r n>1 )) 2 + 1/2 

,pec stub /mT 

°io — J KU, \ 

V e i _ 

+ 

§ (v^TTTT - 1 ) (A/- (v^MT^r"’ 1 ) (v^TT^r"’ 1 )) 2 - 1/2 

(4.72) 

and the more compact representation of the no-stub diffraction coefficient is at once 
evident. 

To assess the importance of the constant &i 0 with respect to Z even (cos 0, cos 0 O ) 
as given in (4.56), we plotted in Figure 4.6 the backscatter echo width patterns 
associated with the three configurations in Figure 4.5(a)-(c) and have compared these 
patterns with that computed by setting = 0. The chosen relative constitutive 
parameters for the left hand side slab are e = 2 and fi = 1.2, and the entire slab is 
of thickness w = .04A (where here A denotes the free space wavelength). We observe 
that the backscatter patterns are, in general, substantially different, underscoring 
the importance of the constant. Although b{l c stub is nearly zero in this case, it will 
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t: T=.04X, e=2-j.0001 , |1=1 .2 Kl 


PEC'* 


\ 

PEC 


Figure 4.6: H z polarization backscatter echo width for a material insert having w = 
.04A, t = 2 — j.000 1 , n = 1.2 modeled with O(w) low contrast GIBCs 
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Figure 4. 
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: H, Polarization bacbcatt® echo width for a material insert having «, = 

l ’ “ ■ ? - 0001 ' ^ = 1.2 modeled with 0(w 3 ) low contrast PIRP, 

*10 = -.4578 +J.2593, 6„ = -.0408 -*0111, L L 
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not be so when w becomes larger. For example, if w is increased to 0.1 A, a 3 rd order 
low contrast GIBC with terms up to 0(w 3 ) is required for an accurate simulation of 
the dielectric. In this case the constants 6 10 , b n and 6 20 are non-zero and as shown 
in Figure 4.7 they play a major role in providing the correct diffracted field by the 
recessed slab of Figure 4.5(a). 

4.3.4 Diffraction by a thin dielectic/ferrite half plane 

Another configuration whose diffraction has been examined in the past is the thin 
dielectric half plane shown in Figure 4.4(c). If a 2 nd order GSTC with terms up to 
O(w) is employed for the simulation of the layer, the resulting diffracted field is given 
by (4.57) and (4.58) in conjunction with the “free space limit” entries in Table 4.1 
and Table 4.2 for the right hand slab. The constant &k> must again be specified for a 
complete determination of the diffraction coefficient. We remark, however, that if 
is arbitrarily set to zero, then for w, = w the sum of (4.57) and (4.58) reduces to the 
diffraction coefficient already derived in [54]. Similar assumptions about the value 
of the constant have also been made in [32] and although this may be acceptable in 
some cases (i.e. for extremely thin layers or layers of certain composition), it was 
already demonstrated above that the constant(s) play an important role and must 
be accurately determined. This is the subject of chapter five. 

4.4 Conclusion 

In conclusion, we have derived a general solution for the diffraction by a discon- 
tinuous distributed sheet representing a multilayered slab discontinuity. The solution 
can be specialized to a wide variety of material junctions and discontinuities by an 
appropriate choice for the polynomial operators U] } and Unfortunately, un- 
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known constants arise in the final solution due to the non-uniqueness of the finite 
GSTC/'GIBC sheet representations. In comparing the derived diffraction coefficient 
to a piev ious result for the material half plane, it was shown that the single unknown 
constant obtained herein was implicitly assumed zero in [o4] and [32]. The impor- 
tance of the constant was examined for the case of three grounded slab geometries 
whose GIBC modeling differed only by the value of a single constant 6 10 . This com- 
parison demonstrated that for very thin metal-dielectric junctions, the constant 6 10 
was approximately zero, whereas for thicker junctions the constant(s) played a more 
crucial role. 



CHAPTER V 


RESOLUTION OF NON-UNIQUENESS 
ASSOCIATED WITH THE GIBC/GSTC 

SOLUTION 


In chapter four, a dual integral equation solution was presented for the diffrac- 
tion by a multilayer material-to-material junction using a GSTC simulation of the 
multilayer slab. As expected, the solution was in terms of unknown constants and it 
was shown that these are dependent on the physical properties of the junction (see 
Figure 5.1). Consequently, an approach for determining the solution constants is to 
enforce tangential field continuity across the junction. This, of course, demands a 
knowledge of the fields internal to the discontinuous slab, which, however, are not 
readily available when a GSTC simulation is employed. The Weiner-Hopf solution 
in conjunction with the GSTC provides only the field external to the slab, and the 
majority of this chapter deals with the determination of the internal field from the 
external one given in chapter four. 

In the following section, a modal representation of the internal field is proposed 
comprised of discrete and continuous spectral components. This representation is 
compatible with that given by Shevchenko [45] whose eigenfunctions are chosen to 
satisfy the continuity boundary conditions across all layer interfaces including the 
air-dielectric interface. Consequently the representation is valid inside and outside 
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Figure 5.1. Sjmmetnc multilayer slab discontinuity illuminated by a plane wave 

the dielectric once the coefficients of the modal representation are determined. This 
is accomplished by recasting the Weiner-Hopf or dual integral equation solution of 
the same problem from chapter four in a form compatible with the proposed modal 
representation, thus permitting the identification of the modal or eigenfunction co- 
efficients. These are, of course, in terms of the unknown constants appearing in 
the Weiner-Hopf solution and the enforcement of field continuity across the junction 
leads to a linear system of equations to be solved for the constants. In the final 
section of the paper, several scattering calculations are presented for a few material- 
to-material junctions which demonstrate the importance of the constants and the 
accuracy of the solution. 



103 



Figure 5.2: Symmetric multilayer slab with irregular termination illuminated by a 
polarized field excitation 

5.1 Modal Decomposition for the Symmetric Slab Problem 

Consider a symmetric slab of total thickness r with an irregular termination to 
its left, as illustrated in Figure 5.2. The slab is herewith assumed to consist of L 
homogeneous layers with the m tfl layer being of thickness r m and having relative 
permittivity and permeability e m and respectively. When this truncated slab is 
subjected to some polarized field excitation, the field to the right of the junction at 
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the point x can be written as a sum of odd- and even-symmetry fields. That i 


is 


F(x,y) = F odd (x,y) + F^ n (x,y) 


( 5 . 1 ) 


where 


E 2 , E z polarization, 

F = l (5.2) 

Z 0 H z , H z polarization, 

F odd (x, y) = -F odd {x,-y) and F' ven (x,y) = F even (x, -y). Following [45], the odd 
and even fields everywhere interior and exterior to the slab may be decomposed into 

discrete and continuous eigenmodes as 

N go N°* d 

F odd (x, y) = £ A ^odd ^ X gof ^ e -jkx\% + £ B odd $ odd (y) e~ jklX ™ d 

m = l m = 1 

ro o 

+ j Q C odd (0) $° dd (A 2 , y ) e~ jkxX d@ (5.3) 

F even (x,y) = £ ,4- en 4' even ((A^) 2 ,y) (if) e"***""" 

to = 1 m = l 

+ C even (0) T even (A 2 , y) e-^d.5 (5.4) 

where /m{A^ d ’ even } < 0 and A = \/T — /3 2 , with the branch of the square root chosen 
so that /m{\/l — T 2 } < 0. In (5.3) and (5.4), \J/ odd - even are referred to as the cross 
section functions corresponding to the continuous modal fields whereas even are 
the corresponding cross section functions for the discrete modal fields associated with 
the surface waves. The cross section function associated with the geometrical optics 
fields is also $«<*<<, even eva j uatec i at X _ ; \go^ w j iere x%[ s a parameter to be determined 
later. As can be observed from (5.3) and (5.4), the cross section functions specify the 
field behavior in the plane normal to the slab, and hence all information pertaining 
to the fields interior to the slab are embedded in these functions. They will be 
chosen to satisfy the orthogonality relations (where u(y) is fi(y) or e(y) for E z or H z 
polarization, respectively) 


*(A 2 ,y)1'(A 2 ,y) 


u {'j) 


dy = 0 for A ^ A 


( 5.51 
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f” *~WJf-y) dy = o (5.6) 

J -> O u(y) 

and thus each discrete eigenmode (y) e~ jkxX and each continuous eigenmode 
'P (A 2 , ij) e~ jkxX must satisfy the wave equation. Additional details pertaining to 
the cross section functions are given in [45]. 


5.1.1 Exterior Cross Section Functions 


To compute the cross section functions in the exterior slab region \y\ > r/2 , we 
recall that in accordance with the slab simulation based on the generalized sheet 
transition conditions (GSTCs), the external fields satisfy the conditions (see chapter 
two) 


l<n (-^) {F m '+} + (-^) {a V F M *} = 0, x>0 (5.7) 

U 2l (-^0 {dyF’"”'*} = 0 . x > 0 . (5 $) 


where F + = F{x,y = r/2), dxF + = £F(x,y = r/2), and dyF + = ^F(x,y) \ y=T / 2 . 
The operators U{ } (— ^5r) are polynomials in —dx 2 /k 2 and their explicit forms for the 
multilayer slab are derived in chapter two as a product of the functions 
<jf,j (u m , —dx 2 / k 2 ) which are completely dependent on the properties of the 

m th layer. In these, the parameter u m is y. m for E z polarization and e m for H z 
polarization, and n m = is the refractive index of the m th layer comprising 

the slab. Because of the orthogonality conditions (5.5) and (5.6), each of the cross 
section functions ^ (A 2 , y) and (y) must satisfy the odd or even GSTC (5.7) or 
(5.8). In view of this 




odd 


( A2 ’ y ) = 17 


sin 


k(\y\-r/2)VT^ 




+ZYi 2 (a 2 ) cos k{\y\ - r/2) Vl-A 2 } 
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= y {l<u (A 2 ) 9 12 (1,1,|j/| - r/2, A 2 ) 
+U X2 (A 2 )«/ 22 (l,l,|y| -r/2,A 2 )} 


^ even (A 2 ,2/) = < U 21 (A 2 ) 


sin k(\y\ - r/2) v^l-A 2 


vT^a 2 

+U 22 (a 2 ) cos fc(|y| - r/2) \/l — A 2 ] } 
= [U 2 \ (A 2 ) < 712 ( 1 , 1 , |t/| - r/2, A 2 ) 

+W 22 (A 2 ) <7 22 (1, 1, |j/j — r/2, A 2 )} 


(5.9) 


(5.10) 


where i/ij (A 2 ) are the same polynomials appearing in (5.7) and(5.8), and q l} represent 
the infinite order form of the layer operators given in chapter two. Once each of 
the modes comprising (5.3) or (5.4) is substituted into (5.7) or (5.8), respectively, the 
differentiation implied by -dx 2 /k 2 reduces to a multiplication by A 2 and the above 
ty odd and 'P et ' en are then readily shown to satisfy the associated GSTC. 

The cross section functions ty odd (A 2 , y) and '$ even (\ 2 ,y) may also be rewritten in 
the form 


ty° dd ^A 2 , y'j = — ^ d ^A 2 ^ e ~2 k (\y\~r/2)yyi-\ i 


f\ 2 \ e J*(M- r / 2 )\/l-A s 

y ^ ' 

(5.11) 

'J' even (A 2 ,y) = ^ ven e -i fe (lv|— r/ 2 )x/l-A J 


+<f et,en (a 2 ) e J fc (lv|-’-/ 2 )v'l-A 2 

(5.12) 


for | y| > r /2, where it can be easily shown that (5.11) and (5.12) satisfy the orthog- 
onality relations (5.5) and (5.6). This representation is customarily employed for the 
surface wave cross section functions. In particular we set 

*£“(») = ; lyl > r/2 (5.13) 

*""(!/) = ; |y| > r/2 (5.14) 
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where A e m ven '^ must now be chosen so that (5.13) and (5.14) satisfy their associated 
GSTC. By substituting (5.13) and (5.14) into (5.7) and (5.8), respectively, we find 
that \^ d ’ even must satisfy the polynomial equations 


7 - wF> 2 U,2 ([a^] 2 ) +W„ ([v") 2 ) = o (5.15) 

1 - (A*ri (K"f) +«21 ([A"'"] ! ) = o. (5.16) 



5.1.2 Interior Cross Section Functions 

We consider now the determination of the cross section functions for the region 
interior to the slab (i.e. in the region |j/| < r/ 2). For simplicity let us first assume 
a single layer slab of thickness r = ‘2TJ, whose upper face is located at y = r t . 
In accordance with the preceeding, the cross section functions associated with the 
external fields are given by 

'i odd (^X 2 ,y) = ^ j< 7 n (^Ui, «!, T x , A 2 j 912(1, 1 , |i/| — ri, A 2 ) 

+?12 («i,«i,T 1 ,A 2 )922 (l,l,|y| -ri,A 2 )}; \y\ > Tj (5.19) 

^ et,en (A 2 .//) = {<721 («i,ici,ri,A a ) <712(1, l,|y| -n,A 2 ) 

+722 (uii «i, iq, A 2 ) <722(1, 1, |y| — Tj, A 2 )} ; |j/| > r, 1 >._' n 
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obtained by setting U tj (A 2 ) = q tJ (u x , k x , t x , A 2 ) in (5.7), (5.9) and (5.10). These are 
orthogonal functions and each must, therefore, satisfy the continuity conditions 


(A 2 , r,-) = <l‘ dJ (a 2 , rf) ,5.21) 

(a 2 , rf) = dy* M (A 2 , rf ) (5.22) 

«*”” (A 2 , rf) = (A 2 , rf) (5.23) 

(A 2 , rf) = Sj,*'*" (A 2 , rf ) (5.24) 

(due to symmetry it is not necessary to enforce similar conditions at y = -t x ). It is 


now straightforward to deduce that possible cross section functions satisfying (5.21) 
- (5.24) are of the form 

'V odd (A 2 , y) = y?i 2 (ui,Ki,|y|,A 2 ) (5.25) 

y even ( A2 ,y) = ?22(«i,«!,|y|,A 2 ) (5.26) 


for |y| < Ti. Also, in view of (5.21) - (5.24), the cross section functions for the surface 
wave modes remain as given in (5.17) - (5.18), provided (5.25) and (5.26) are used 
in place of 

As a specific example, let us consider a low contrast 0(t x ) representation of a 
single layer. From chapter four, the truncated GSTC operators IA %] ( — dx 2 /k 2 ) are 
given to 0(t x ) by 


U\\ 

U\2 

Ui\ 


dx 2 ' 

' 

dx 2 ' 

k 2 , 


= U 


n 


dx 2 ' 

k 2 


= 1 


= jku x r x 


dx 2 \ _ jkTyK 2 


k 2 


Ul 


+ 


p-dx 2 

kxii 


and the corresponding q tJ (A 2 ) polynomials become 


(5.27) 

(5.28) 

(5.29) 


9 xi («i, «i, r x , A 2 ) = 922 («i,«i,Ti,A 2 ) = 1 


5.30) 
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and 


<712 (ui,Ki,ri, A 2 ) = jkuiTi 

?21 (ui, /Cl, Ti. A 2 ) = — A 2 ) 


(5.31) 

(5.32) 


912 ^«1,«1, |j/|, A 2 ) = jAruxlyl 

922 (ui,«l, M, A 2 ) = 1. 

When these are substituted into (5.19), (5.20), (5.25) and (5.26) we obtain 


(5.33) 

(5.34) 


<H odd (\ 2 ,y) = — 

v ' y 


. sinffc(|y| — Tj )\/ 1— A 2 

J ‘ /nr — 


+;fcu 1 r 1 cos fc(|y| - Ty) \/l-A 2 ] ; |y| > ( 5 - 3 °. 


jkuilyl; 


\y\ < Ti 


^> sven (X \y) = 


-fcrj ( 2 \ 2 \ sin(fc (ly|-n)Vl-A 2 

Ul (*1 “ A ) 


(5.36) 


+ cos [k{\y\ - Ti) v/l-A 2 ; |y| > n 

{ 1; |y| < T X 

as the cross section functions for a single layer slab modelled with an 0{r x ) GSTC. 

For the general case of a multilayer slab, it is necessary that each of the internal 
cross sections functions satisfy the continuity conditions at all layer interfaces com* 
prising the slab. From chapter two, we obtain that the boundary conditions at the 
interface between the L — I th and L th layer are 


V L ' 


L-l 


dx r 
' k 2 


{u lF““(x, !, = !,£_,)} 


V, 


x - 1 
21 


■jr) {«L/""”(x, !, = !,£_.)} 


]_^nL — 1 

k 


+ tV 22 


dx 2 

"W 


dx 2 

k 2 


{dyF odd (*,y = y£_x)}=G. 37) 


{dyF™ (x,y = y£ l )}=(W 8) 


where xji is the y coordinate associated with the top surface of the L th layer and the 
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operators V fj 1 (— are given by 


r Tfr 1 

f ) 

{ k* ) 

i Pn 1 

( dx 2 \ -> 

V * 2 J 

L-l 

=n 

m = l 

^ m i 7 ^2 ) 7 1 2 ( ^ tti ? ^ tti 7 7 ^ ) 

. -p.fr 1 1 


1 

f i)r 2 \ 

^ ?21 (^7717 ^ 771 7 ^2 ) 922(^7717 ^ 771 7 ^ 771 7 "p“ ) J 


Possible interior cross section functions satisfying these conditions are 

*“"(A \y) = ld{ c g( A 2) ?1J („ LiSt> | ! ,|_ ! , t _ ltA 2) 

+-C* 2 (a 2 ) 922 (^l, /cl, |y| — j/L-i) A 2 ) j vl-\ < |y| < yi5.40) 
'P eve "(A 2 ,y) = {x?rrM«i 2 («i.*L,|y|-yL-i,A 2 ) 

+X>aT ( A2 ) 9m(“£> /Cl, |y| - y^-i, A 2 )} y L _j < |y| < j(p.41) 

in which £) even ' odd (A 2 ) are to be determined by demanding that the interior and 
exterior cross section functions are continuous at y = r/2(= Ul)- Setting y = r/2 in 
the expressions for the external cross section functions (5.9) and (5.10), we obtain 

ty odd (A 2 , y = r/2) = U x2 (\ 2 ) (5.42) 

* euen (A 2 ,y = r/2) = W 22 (A 2 ) (5.43) 

and when these are equated to (5.40) and (5.41) we find 

4°" (A 2 , d = 7 {^U 1 (a 2 ) 91, K. .t, Ivl - A 2 ) 

+ f’lf 1 (A 2 ) 9„(“I.1 «i> lyl - 91 - 1 . A 2 )} (5.44) 

*""(A 2 ,y) = V$T l (A 2 ) k l , |i/| - y L . u A 2 ) 

d"^22 (A ) 922 ( u Li /C Li |y| - J/L-l, A ) (o. 45 ) 

for y L > \y\ > y L . x . 


To derive the cross section functions for the other layers, the above procedure 
may be repeated in a recursive manner until all layers are accounted for. Doing so. 


Ill 


we find that a complete representation of the cross section function everywhere is 

r (A 2 ) 9 j 2 (1, 1, \y\ - r/2, A 2 ) 

+Wj 2 (A 2 ) <722 ( 1' 1' ll/l - t/2,A 2 ); \y\ > r/2 

« 0 "(A‘ 2 f j/) = y ^{r 1 (A 2 )?i 2 («/,«/,l</|-'//-i,A 2 ) (5.46) 

+v[i 1 (A 2 ) 922(«/i */, li/l - yi-i, A 2 ); yi > |y| > yi - 1 
. 9w(«i,«i» I 2 /I 1 A 2 ); \y\<yi 

' U 21 (A 2 ) 912 ( 1 . 1, lyl - 7-/2. A 2 ) 

+Z ./22 (A 2 ) 922 ( 1 , 1, jy| — r/2, A 2 ); |y| > t /- 

* even (\ \y) = I ^r 1 (A 2 ) 9 i 2 (tt/,K ; ,|y|-y ; _i,A 2 ) (5.47) 

+^22 1 (A 2 ) 922 ( 77 /, Ki, \y\ - yz- 1 , A 2 ); y { > |y| > i 
. 922 ( 771 , «i,|y|, A 2 ); \y\<yi 

when these are used in (5.3) and (5.4) in conjunction with (5.17) and (5.18) we have 
a complete field representation for x > 0. 

5.2 Recasting of the Dual Integral Equation Solution for a 
Material Discontinuity 

The expansions (5.3) and (5.4) can be used to represent the fields interior and 
exterior to the discontinuous slab shown in Figure 5.1. For x > 0, the material 
parameters used in the definition of the cross section functions (5.46) - (5.47) must 
then be associated with the right hand portion of the slab. Similarly, for x < 0, the 
material parameters in (5.46) - (5.47) must be those of the left hand portion of the 
slab. 

The diffraction by the slab discontinuity shown in Figure 5.1 was the subject of 
chapter four where a complete expression for the scattered field was given by em- 
ploying a GSTC simulation of the slab. However, owing to the non-uniqueness of the 
GSTC, the resulting diffraction coefficient was in terms of unknown constants whose 
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determination was shown to require additional constraints. One such constraint is 
the enforcement of a boundary condition demanding field continuity at the material 
junction, but this requires a knowledge of the slab interior fields and the solution 
given in chapter four pertains only to external fields. However, as shown in the 
preceeding section, the expansion (5.3) and (5.4) is valid everywhere when used in 
conjunction with the cross section functions given in (5.46)-(5.47). Moreover, since 
the expansion coefficients remain unchanged for the exterior and interior fields, once 
determined, the representation m.3) and (5.4) can be used to find the field every- 
where. Since the exterior field associated with the slab discontinuity in Figure 5.1 
has already been given in chapter four, it can be used to identify the expansion co- 
efficients. This requires that the solution in chapter four is first recast into a form 
compatible with that in (5.3) and (5.4), making possible the identification of the un- 
known constants which can then be determined by enforcing field continuity across 
the junction. In the following, upon stating the exterior solution we then proceed 
with the identification of the expansion coefficients. 

In chapter four, the discontinuous multilayer symmetric slab shown in Figure 5.1 
was simulated by the GSTCs 

u'u {C“ JA+ } + ] f<n (-^r) { dyF «“■+} =0, -co < x < ((5.48) 

+ =0, -co < x < ((5.49) 

«11 (-^-r) T""'*} + \ul \ =0, 0 < x < CO (5.50) 

K (- jr) {m ,n ' + } + ji<i 2 (ypr) {dj/F*'*" + } = o, o < x < co. ( 5 . 51 ) 

with the superscripts 1 and 2 again denoting the material to the left and right of the 
junction, respectively. The total field is the sum of the even and odd components 
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and can be represented as 


F(x.y) = F odd (x, y) 4- F even (x, y) 


where 


Fine (a:, y) + F,(x,y) + F ref i(x,y ) y > r/2 
Ftran {x,y) + F s (x,y) y < -r/2 

F inc (x,y ) = c J*(*C 0 .*.+y«n* o ) 


'5.52) 


'5.53) 


is the incident plane wave field. F re fi(x,y ) denotes the reflected field, which from 
(4.5) and (4.7) is gi’-en by 

fi jkr sin 0 o 


Frejl V) 


R\ vcn + R\ 


odd 


0 jk(x cos <po-y sin <po) 


(5.54) 


with R\ dd and FL~^ en given by (4.9) and (4.10), respectively. Similarly, F tran , (x, y) is 
the transmitted field, which from (4.6) and (4.8) is given as 


Ftran (•£> 2/) 


0 jkr sin <t> 0 


j^even ^jk{xcos <p 0 + y sin <£ 0 ) 


0.00 


Finally, F s (x. y) is the field scattered by the discontinuity and upon employing the 
dual integral equation method in conjunction with the GSTC (5.48) - (5.51), we find 
(from chapter four) 

F, (x, y) = F° s dd (x, y) + F' ven (x, y) (5.56) 


W 


here 


T^odd { \ V f°° J 

F • (X ’ y) = M-/-oo2 




1 2/| J- oo 2 tt A + A 0 Glt d (A) Q\ dd (A 0 ) G? + d (A) (A 0 ) 


^ odd 1 ^ odd 1 tn 

Z odd (~^o) + Z E «mn(A + A 0 ) m (AA 0 ) n 

m=l n=Q 


e -]k\y\s/i=>? e -jkx\ d ^ 


F, 


r n (x,y) = r - 

J — OO l 


j v/i-A^vT^ 


1 — A 2 


(5.57) 


2tt A + A q Q{ v J n (A) GV- - n ( A o) GZl en (A) GIT (A 0 ) 


1 J Veven — 1 — TTl 

z m n(-AA 0 )+E E ^„(A + A 0 ) m (AA 0 ) n 

m— 1 n=0 


-r=== 1 -0-.5S) 

v/l-A 2 
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In these expressions, X 0 

= cos do, and 


odd 

= int { 1/2 (Nl dd + Nh d + l) } 

(5.59) 

Y even 

= int { 1/2 [Nl ven + N 2 ven + l) } 

(5.60) 

£? + d {\)G°d d (\) 

= Q° x dd (A 2 ) = U l n (A 2 ) + y/mu\ 2 (A 2 ) 

(5.61) 

QtT( A)Grr(A) 

= Gr n (a 2 ) = U\ x (A 2 ) + y/i=XU\ 2 (A 2 ) 

(5.62) 

$3? (A) as? (A) 

= Qf d (a 2 ) = u 2 n (a 2 ) + (A 2 ) 

(5.63) 

g\T (A) g\ v _ tn (A) 

= QT n (a 2 ) = u\ x (A 2 ) + '/T ~u 2 2 (a 2 ) 

(5.64) 

Z 0 dd (— AA 0 ) 

= [^ii ( AA 0 ) 4./ 12 (— AA 0 ) — U\ 2 {— AA 0 ) U xl (— AA 0 )j 

(5.65) 

Zeven ( AA 0 ) 

= [^21 ( A A o) 7^22 ( — AA 0 ) — U\ 2 (~ AA 0 ) U\ x (— AA„)J 

(5.66) 

V 1 

A odd 

= 0(A) of Q\ dd (A 2 ) 

(5.67) 

V 2 

^ odd 

= 0(A) of g° dd (a 2 ) 

(5.68) 

N l 

even 

= 0(A) of C/* v * n (A 2 ) 

(5.69) 

.V 2 

even 

= 0(A) of gr n (a 2 ) 

(5.70) 

As seen, the scattered field expressions are in terms of unknown constants a mn 


and b mn and to determine these via the procedure outlined above, we must first 
rewrite F(x,y ) in a form compatible with ( 5 . 3 ) and ( 5 . 4 ). To do this we need to 
identify from (5.52) to (5.58) the discrete and continuous spectral components. The 
discrete portion of the spectrum is, of course, comprised of the geometrical optics 
and the surface wave fields. These can be identified by detouring the integration path 
in (5.57) and (5.58) as shown in Figure 5.3. In particular, for x < 0 the integration 
path may be deformed to one over the branch cut in the upper half of the A plane, 
capturing any surface wave poles attributed to the zeros of Q°t d ( A) and Q\ v J n (A). 
Similarly, for x > 0, the integration path may be deformed to one over the branch 
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cut in the lower half of the A plane causing the capture of the geometrical optics pole 
at A = — A c in addition to any surface wave poles attributed to the zeros of (A) 
and (A). 

Through the above deformation of the integration paths in (5.57) and (5.58) we 
obtain 


F odd (x,y) = 


F even (x,y) = { 

where for \y\ > r/2 


(x, y) + F}** (x.y) + Fj,f/ x < 0 

y) + F;* dd (x.y) + (x, y) x > 0 

( Fir’" (x, y) + (X, y) + F' u ’/r (X, y) X < 0 

F‘r"‘ (x, y) + F^’ d (x,y) + Flrr (x, y) x > 0 


;s.7i) 


77 1 .odd 
*30 

(*,y) 

r?2.odd 
^ jo 

(x<y) 

rl ( et'en 

**30 

{x-y) 

7^2. even 
^ go 

{x.y) 


j^odd 


Deven 

1^2 

ir\,odd 

t su.- 

(x,y) 


y e jkxCO& 4 >o 

^ | e j^|v|sinqi>o ftodd e - j k(\y\-r) sin <t>oj 

j kx cos <t> 0 

& _ gjk\y \ sin do _|_ f£°dd e -jk(\y\-r) sin do j 


li/l 


ojkx COS <t>o 


| e jfc|y|sin0o A:(|y| -r) sin do j 


0 jkx cos do 


{■ 


j/c|y|sindo jfc(|y|-r) sind< 


■} 


sin <t> 0 U \2 (cos 2 0 O ) — U\ x (cos 2 4 > 0 ) 
sin <pM ?2 ( cos2 0o) + U\ x (cos 2 4 > 0 ) 
sin (cos 2 <f> 0 ) — U\ x (cos 2 <f> 0 ) 
sin 4 > 0 U 2 2 2 (cos 2 4 > 0 ) + Ul x (cos 2 <j> 0 ) ' 


(5.73) 

(5.74) 

(5.75) 

(5.76) 

(5.77) 

(5.78) 


Nl 


— sin <p 0 e 


jl T -jkxXj odd e -j*|y|0-(A‘- 


|y| l=l \j’° dd + A 0 £|dd ( Aq ) £o£d ( Ao ) 


agfl d (A) 


aA 


J A=A: 


1 ^ odd ^ — 771 

z^{-K\'r di )+ Y. E »™(^° , “+a„)"(a;-""a,,)” 

m = l n=0 


(5.79) 


iVf, 


• i J "T* f sin 4>o + \J 1 ~ ( A, ,<>dd jl \2, odd - * | j A / \ 2, odd v2 

FWhr.i-i-r^ 111 ^ c 1 ;e ' e_j |ylv/ (t ’ 

V J ’ 'J) ~ | | > 2,o 


\y\ l=l \*' 0dd + A 0 £o£d ( A 2 -^) £o£«* (Ao) £odd ( Ao ) ££& 


dg° dd (\) 

d\ 


\=\]'° dd 


^ odd ~ 1 N odd 1 Fn 

z M (- E E a mn (x] iM + K) m {\r dd x „)” 

m = 1 n=0 


5. SO) 
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J T i sin 


in 0 O +^1 - (A, 1 )> ) _ jkxX i ...» _jJt lylv /r^(Aj 


1 ,eti«n ^2 


^ A,— + A 0 (a}-—) £|«« (A 0 ) a?" 8 " (A 0 ) - 

A 1 e u « n 1 ^<vert“l — TTl 


*3£?I <n (A) 


a a 


> i l,ei'<n 

A — 'i 


5.S1) 


40™ (-A.A, 1 '"") + £ ' £ 6 w.(A, , -*~“ + A„)"(A, I -*"*A„) n 

m=l n = 0 

^(.in^o + v/l-(*? , " ,B ) 2 ) -jkx\]- nm e -jfc|y| N /l-(A*-« t '*")a 


y2,«t/en , , J 

(*.»)=£ ~ sm0 ° e 


a^T n (A) 

9A 


=x x i ' even + A ° G{ v j n (A 2,et,en ) g\ v j n (A 0 ) g^ v + en (A 0 ) 

Ae ygrt ~ 1 'A« v en — 1 — m , 

4™, (-A„A?''~) + E E 6m„ (A?-’”'" + A.) (A?-'-"A„j 


A= A: 


m = l n=0 


5 .82) 


771 , 0 ^^ / \ y f j 

F d '" (x ’ y)= 


• ,yT^vT^ 6 ;W2 ( v /h+vwi) 


\y\ H 2 tt a + A 0 (A) a?? (A.) 0 2 T (A) (7 S? (Ao) 


Aodd 1 ^odd 1 T7T 

^(-AA 0 ) + £ XI «mn(A + AX(AA 0 ) n 

m = l n=0 


-jk\y\Vi^ e -jkxX dX 

— (5.83) 


f 2 .odd, \ JJ_ f j \A A 2 

dlf} {x,y) ~ \y\ Jct r 2 x A + A 0 g\ 


jkr/2 (^y/ l-Aj+x/l-A 2 ^ 


N odd “ 1 ^cdd - 1 _r n 

Z«w(-AA 0 ) + £ £ «mn(A + A 0 ) m (AA 0 ) n 

m=l n = 0 


odd (\)g° 1 d J(x 0 )g^ d (x)gf + d (\ 0 ) 

e -jk\yUi=F e -]kx\d\ 




5.84) 


n? s v r^y)= J c 


^/l-A^x/I^A 7 e ifcr/2(v^+vOT) 


c£ 27T A + A 0 gf_ d (X)gfj (A 0 ) £ 2 °f (A) g° 2 i d (x 0 ) 


A r «t ;en~*l v « n ~ 1 ~ ^ 

^, e n(-AA 0 )+£ £ 6 mn (A + A 0 ) m (AA o r 

m = l n=0 


,-jk |y|\/l-A 2 e -jfcxA^^ 




65-85) 


F 2 d r;r^y)= L 




cl 2t A + A 0 ar (A) (A a ) a 2 T (A) GaT (A„) 


■^evcn“*l N«t/en ~~ 1 — Tr l 

Z e „ en ( — AA 0 ) + £ £ 6 mn (A + A 0 ) m (AA 0 ) n 

77i = 1 n = 0 


e -]k\y\s/T3J -]kx\J\ 

45-36) 


n/I^A* 


In these, the components F go , F sw , F d i / / denote the geometrical optics, surface wave, 
and branch cut (or diffraction) contributions to the total fields. Also N]^ is the 
number of captured surface waves (i.e., those with normalized propagation constants 
A, 1 ’ 2 such that Im {A; 1 } > 0 or Im {A 2 } < 0). 
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To identify the expansion coefficients A m , B m , and C(A) appearing in (5.3) and 
(5.4), (5.73) - (5.86) must first be rearranged. By making use of (4.9) - (4.10), (5.61) 
- (5.64), and (5.77) - (5.78), the geometrical optics terms can be rewritten as 


F. 


1 .odd / 
go 


(x.y) = 


F 


2 ,odd 
3 o 


x.y) = 


Fh' ven {x,y) = 


30 


p2 ,even 


30 


(*.y) = 


sin d> 0 6 jkxcos <t>o e jkr ^ 2 sm< ^° 


\v\ 

6t dd ( a 2 ) 


-rU\ 2 

(A 2 ) cos ( k sin 

bo [|y| 

y sin 6 

fi jkx cos 4>o~jkr/2 sin^o 
qC c 

li/I 

Q°2 dd ( A 2 ) 


+£Vl2 

cos (k sin 

<t>o[\y\ 

y sin <& 

e jkx cos <b 0 g jkr/2 sin <i> Q 


ar n (A 2 ) 



^A 2 j cos ( k sin 

My} 

y sin 0 

#jkx cos 4> 0 ~jkr/ 2 sin <t> 0 
O 




Ul (a 2 ) 


Ul (a 2 ) 


C?“( K) 

+^22 (A 2 ) cos (fcsin^ [|y| - r/2])] 


2 \ j sin ( k sin d 0 [|y| - T l-\) 


sin <p 0 

(5.87) 

j sin (k sin <t> 0 [|y| 

- r/2}) 

sin d 0 

(5.88) 

j sin (ksindo [|i/| 

-r/2}) 

sin do 

(5.89) 

j sin (k sin do [|j/| 

-r/2}) 


sin d 0 


(5.90) 


The bracketed terms in (5.87) - (5.90) are now readily recognized as the cross section 
functions given by (5.46) and (5.47) once \ 0 is set to cos do- Thus, from (5.3) and 
(5.4) the geometrical optics fields may be expressed as 


F g l ? dd (x,y) = A\ d \\ 0 )V\ dd (\l, y )e> kx ™*° 
Fg 0 ° dd (x,y) = A° 2 dd (\ 0 )V 0 2 dd (\ 2 0 ,y)e> kxcoa *° 

F]r n (*, y) = A\ ven (A 0 ) (A 2 , y) e’ kx ™*' 

F 2 g r n (*, y) = A? en (A 0 ) *r n (A*, y) e> kx ™*' 


where the A expansion coefficients are identified as 

sin 


A^iX) = 
A? d (A.) = 
.-ir'VA.i = 


er“(A 2 > 

sind> 0 e-* fcT/2sin ^ 

WlK) 

sin ^ oC i^/2«in^o 

^r en (A 2 ) 


(5.91) 

(5.92) 

(5.93) 

(5.94) 


(5.95) 

(5.96) 

(5.97) 
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A? en (A„) = 


sin (f) 0 e^ kT l 2s ' n ' t>0 

Gr n iK) 


(5.9S) 


Also. #f“* ei ' en are the same as those given in (5.46)-(5.47) but are associated with 
material parameters to the left of the slab and likewise \jr° dd ’ even are associated with 
the material parameters to the right of the slab. 

To identify the B expansion coefficients we observe that the surface wave terms 
are already of the form appearing in (5.3) and (5.4) and can be readily rewritten as 

vrl ,odd 


F]: dd {x,y) = £ B) Md {K) {Z 0dd (-\ 0 \]'° dd ) 


1=1 


^ odd 1 odd 1 171 . m v „ 

+ E y. +A ») ( A ? A °) 

m=l n=0 


/ f ,odd 


■ i .odd 

~ jkl \ (5.99) 


,y\ c 


J V 3UI 

F%“{x,y) = Y \Zodd (—Ktf ,odd ) 


1 = 1 


+ (A?'»" + A„) m (a?'“"A„)" *?“ ([A?'”"] 2 ,!,) 


m = 1 n = 0 


3 -jkx\] 


2. odd 
l 


5.100) 


»d,et<n 

3 


; j ven (x,y) = £ Bl even ( X 0 )[ Z 0 dd (- X 0 \ l r en ) 


1=1 

A r even — 1 .V tflfen — 1 — rn 




+ Y E bmn (A, 1,eU * n +A 0 ) m (A ( 1 ’* Ve ' , A 0 ) 

m = l n=0 

^2, even 

C,y"”(x,y) = £ (A.) \Z M (-A„A?'"“) 

/=1 

+ £ £ 6m, (A?~“+ A c ) ” (A?-~"A„) * r"( [A?' 

m = 1 n=0 


e„ e nf r 1 l.ewe»l 2 y j g-J^^foi) 


,y e 


-jkx 


\J,<uen 

^5.102) 


with 


B]'° dd (\ 0 ) 


Bl odd { K) 


sin 4>o 


c jkT/2sin<t> 0 


\]' odd + A c 
— sin (j> 0 


a^i d (A) 
aA 


J A = A 


Sit i (A.) (a}' 0 ") esf( A„) 

gjkr / 2 sin 4> 0 


(5.103) 


A 


2, odd 

l 


+ x o g° d J (a l° dd ) g° x t d (A 0 ) 


dC%i d (X) 

a\ 


i A=A*-°" 


GF + d (K) 


(5.104) 
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B}' even (A„) = 


— sin o 0 

A p en T Aj 


gjkr/2 sin <i 0 


A . A ,.». n <*.> sir (a;—) ggr ( A, 


5 ; 


Aj = 


~ sin o 0 

Af— + A c 


gjkr /2sin<i> 0 


5.105 i 


cpr (A, 2 '~")ef».-(A„) 


gg^lLll 

d\ 


A=l . (a.) 


' 5 . 106 ) 


It 


remains to identify the C expansion coefficients and to do this, it is necessary 
t° employ the transformation J = vT^. Doing so in (5.S3) permits ns to rewrite 


r~ l ,cdd 
f dt/f as 


Flffix.y) = 2- i ~ 3 \0 Z K 

1 .V J — x> 2- . /I 22 ^ /t 77 


^ 2 * 7r- j 2 (ym* + a„) 

; ; 1/2 ('v/ 1-A o+' 3 ) e |<3 e 


^ ( n/T ~ jj + * yf) g°i- ( A °) a 2 °f (yi^F) G$? (A 0 ) 

__ / / \ ^ r odd - 1 Wpdd — 1 — rn 

•“ E (A 0+v ^) m ( Ao yrrjij” 


'<(3.107) 


"here the branch of the square root is chosen so that Im(JT=F) > 0 and 6 is a 
vanishingly small positive number. The integral may now be split into its positive 
and negative portions and then recombined to obtain 

_l£i Hz3 




Tl / 7 ; — ua u V v ' 5 /x o \ 

W 0 2e V Tr3i( v r T r ^+A 0 )ef" ( A o) e!"(vTr^) £ , |?(Ao) 

e-^(|y|-r/2)/3 eiMW-r/aw ) 


. (v^F + <5) (VTFF - *) 


d0 


(5.108) 


where 


z M (sj r^.A.) = z^(-A„yr^) 

•^odd - 1 ^odd — 1 — m 

+ 5 £ Gmn ( A ° + (^a/^f) 71 (5.109) 
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By invoking (5.61), vve can show that 


:-jk(\y\-r/2)3 


ojk{\u\-T/7)3 


-23G°j d (yT^J 7 ) 


Q\ 7- d + <5) QfJ ( - S) [Uh ( l - 3 2 )} 2 - 3 2 [£A\ ( l - ^)j : 

7 — (l - 3 2 ) + cos (fc£ [|y| - t/ 2]) Wj 2 (l - ,^ 2 ) (5 


5.110) 

where the term in square brackets is identified as ^^(l — 3 2 ,y) as given in (5.46) 
for \y\ > r/2. Substituting (5.110) into (5.108) finally yields 


F I .odd / \ 

diff (x,y) 


C odd ( J) 


•'odd 


i-3*,\ 0 )'H° dd (i-3 2 ,y)e 


-jk XS /i-pi d3 


(5.111) 


with the expansion coefficient C° dd (3) given by 


C° dd (3) = ~ 


3\^ 


* vT^l 1 (x/T^T 1 + A 0 ) 

e :kr,2^ 0 godd (yr^) 


Qft (A«) (A 0 ) {[Uh (l - 3 2 )) 2 - 3 2 [^ 2 (l - 

In a parallel manner we obtain 


*frll2) 


r-r2,OcW / \ 

^t// (x,y) = 


' C° 2 dd (3) Zodd (\A^, A 0 ) n dd (i - 3 2 , y) e-^ x ^d3 

(5.113) 

Fdi}7 n (X, y) = / 0 °° cr en (/?) 2 even (yr^ 2 , Ao) ^r n (l - /? 2 , 2/) e~ jkx ^d3 

(5.114) 

Fd>f V } en (x, y) = ^°° cr n (/?) 2 even A 0 ) ^r n (l - £ 2 , y) 

(5.115) 

where 


l-^ 2 ,A,j = Z' Ven (-\ 0 Jl^3^ 

^«fen“l ^Vgven — 1 “ Wl / . \ m /■ \ /t 

+ E E 6mn(A.+ yi-^) A o yi-^)5.116) 


m=l n=0 
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and the spectral expansion coefficients are given by 


Cf d (3) 


Fy/l-Vo 


* v/T^F + A ( 


0 jkr/2Wl^l r odd 


GfJ (x/T^F) Q\- (K) G& d (A„) {[U 2 n (1 - /3 2 )] 2 - 3 2 [Ul 2 (l - J 2 )]' ? y ' 

„ = l NE*L 

T v'T^(v/T^^A 0 ) 

e jkTl2^igodd ( 

(yr^) $§" (X 0 ) g*t d (A 0 ) {[W 2 \ (1 - ¥)\ 2 - p [U& (i - 3^f lb] 


,3\ l-M 


-3 2 + A 0 


e jkr/2yU^igodd 


C 2 eLen (3) = -i 


QfJ (s/i=F) G{t d (A 0 ) (A 0 ) {[U 2 y (1 - fl 2 )] 2 - 0 2 [U&(1 - 0*)YY 

From (5.3) and (5.4), the modal terms (5.91) - (5.94), (5.99)-(5.102), (5.111), and 
(5.113) - (5.115) provide a field representation which is valid everywhere. Since the 
odd and even fields are decoupled, two independent representations are obtained for 
each of these fields. 


5.3 Determination of the Constants 


To determine the constants a mn and 6 mn , we may now enforce the tangential field 
continuity conditions 

F(x = 0“,i/) = F (x = 0 + ,y) ; |j/| < r/2 (5.120) 

^axF(x, #)„ o- = z7M dxF(i ' ! ')*= 0 * ; |!,i<T/2 (- 5 - 121 ) 

with 

(m 1.2 (y) £z-pol 

«u(y)= < (5.122) 

Ui, 2 (y) H z - pol 
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and the subscripts 1 and 2 denoting quantities attributed to the left and right side 
of the slab. Substituting (5.71) - (5.72) into (5.120) and (5.121) , we obtain 


= F^ dd (x = 0 + ,y) 


F}„ odd (x,y) + F' d ;f f d (x,y)} x=Q _ 
l odd {x,y) + F]; 0 // (x,y)} £=0+ 
X di // IT 


■dx \F^ odd {x,y) + 

^2, odd ( \1 

d.ff {x,y)\ 


My) 

C" (x = 0-,y) + F,\; L,n (x = 0~,y) + 

0 + ,s<) + F’r" (* = 0*.») + FJ5T t 1 


J?2,tven t 
go \ 


— t-t5x 

Ul (y) 

pi, even / 

L 30 

— ! — dx 

F?f ven (. 

U2 {y) 

go v 

for all a mn and b 


0 ,y) 

0 + ,y) 

F? w :ven (*» v) + F d,T/ n ( x -y). 


: 0 - 
x=0 + 


(5. 123) 


(5.124) 


(5.125) 


(5.126) 


to be solved for all a mn and b mn . In particular, for an odd GSTC of 0{N{ 


dd \ 


to 


tiie left and of 0(iV| dd ) to the right of the discontinuity, the number of a mn to be 
determined is equal to 

(\7 i\ f ( N ldd+^dd)(^ldd+^dd~ 2 ) . \ri | v 2 * c - cn 

-'odd [ *odd - i J a S + iV oci£i is even 

v > ( 0.12 () 

JViu + M is °dd 


■Y„ = 


( , Vdd'*~ j ’''odd) 1 . 

8 ’ 

To determine all a constants, (5.123) and/or (5.124) must then be enforced or sam- 
pled at a minimum of N a points across |y| < r/2 and 0 < <t> 0 < 7r. Similarly for an 
even GSTC of 0{N{ ven ) to the left and of 0(A r | ven ) to the right of the discontinuity, 


Neven i^^’even l) 

V * = 2 


f ( Nj vm + *%,,„ ) ( NU m - 2 ) . v 2 ; 

a i iV even ' iy euen 1 


= < 


is even 


(Njven + Njven ) . 

8 ’ 


+ A ’elen » odd 

(5.128) 


and thus, the b constants can be determined by enforcing (5.125) and/or (5.126) at 
a minimum of A5, points. 

Substituting for the fields in (5.123) and (5.124) as given in the previous section. 
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Figure o -i. indexing scheme for constants. 


we obtain the equations 

V^(K,y) 
Vg d F ( K,y) 

where a p = a m(p ) n ( p) with 


iV 

Yl a pQF d (m(p) 


P= 1 


5 ^ (p) ? ^0? 2/ ) 


N a 

Y1 a pQdxF ( m (?) .n(p), A 0 , y) 


P=1 


(5.129) 

(5.130) 


P 

m (p) 
n (p) 


[n + m — 1) (m + n) 


4- m 


P — — 7n< 


i,-J v' 1 + S(p-l)-l l ja f v /l + 8(p-l) + l 


7n£ 


y/l+8(p-l) + l j 


- m(p) 


(5.131) 

(5.132) 

(5.133) 


which are in accordance with the ordering of the a mn constants as the order of the 
GSTC is increased (see Figure 5.4). Also, 


ATl,0dd 

JV JW 

+ £ B]' oid (\ 0 ) Z 0 u (— AA/ -0 ^) 'b° dd ([A ( 1,wM ] 2 , y) 

N 2 , w ° dd 
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Provided X°j d are even, the integrals in (5.135) and (5.137) converge (see Table 5.1) 
and (5.129) - (5.130) can be solved for the constants a p by matching (5.129) at A'„i 
points and (5.130) at N a2 points such that N al + N a2 > N a . 

To solve" for the b p constants we substitute for the fields in (5.125) and/or (5.126). 


giving 

‘ V 6 

Vr'(K,y) = EWr"("i(p),n(p).A„,j) (5.13S) 

P=1 

Nb 

v dx F n (K,y) = (5.139) 

p=i 

with 


t /even 
V F 


(x 0 , y ) = Ar n (x 0 )n ven (xiy)-Ar n (K)n ven (^y) 

y 1 , e 1 1 e n 

+ 22 Bf’ even (A 0 ) Z even (— AA^’ even ) f [A / 1 ’ et ' 


' y 


VaTFi^y) = 


- 22 B f' even (■ K) Z even (~\\? even ) ( [Af ,et '* n ] 2 , y 

+ f~CV m WZ mn (\- 3 \y)di 3 

~ £° Cr" (?) (-/-PK'j (l - 0\y) d.3 

‘ P"’"'' - u 2 (y) 2 V K ’ y > 


45.140) 


«i (y) 


N\ 


1 ,ev«n 


- E 

i=i 


Aft 


+ E 


\w**Bl' even (\ 0 ) 

ui (y) 

\ 2, even r>2,et>en / \ \ 

\ I'M 


Zeven (~XX]' eVen ) 


even / f \ 1 .even] ^ 


/=1 


u 2 (y) 


Z eV en (- AA^) ( [X^f ’ , 


■ I — Tw - (-^.) (1 - /J 2 . S’) d, 3 

[<* y/l — /3 2 C* ven (B) / / \ , 

Jo (-^) (‘ - ^ ») " 

(5.141) 


+ 


and 


QT n (K.y) 
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S° x dd or 

'N° dd or 

Junction 

Branch cut integral 

Order of branch 

yeicn 

1 

\'eren 
- N 2 

Field 

(odd or even) 

cut integrand 

even 

even 

E z or H z 

Fhss( x = 

\p\- 3 




Fl n (x= 0+,jr) 

i/r 3 

odd 

odd 

E z or H z 

E} i t f (x - 0 _ . y) 

m~ 2 




F] lff {x = 0 + ,y) \ 

\3 \ - 2 

even 

odd 

E z or H z 

1 

O 

II 

H 

S? 

\I3\~ S/2 




Flffix = 0 + ,y) 

|,d|" 3/2 

odd 

even 

E z or H z 

F}iss ( x = o i y) 

\t>\- w 




Fjiss ( x = °Vy) 

l/3|- J/2 

even 

even 

H y or E y 

F d,ss( x = °~’») 

m- 2 




F] iff {x = 0 + , y) 

\@\~ 2 

odd 

odd 

Hy Or Ey 
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Table 5.1: Asymptotic behavior of integrand for the functions Fuss ( x = 
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(5.143) 


As before 6 P = 6 m ( p ) n ( p ) with p, m(p), and n(p) given in (5 . 1 31 )- (5 . 1 33) . Again, 
provided jV$ en are even, the integrals in (5.141) and (5.143) converge and (5.138) 
- (5.139) can be solved for the constants b p by matching (5.138) at N bi points and 
(5.139) at N b2 points such that N bx + N b2 > N b . Results based on the solution of 
(5.129) - (5.130) and (5.138) - (5.139) are considered next. 

5,4 Validation of the Solution 


In this section we address the validation of the GSTC solution. In particular, 
several diffraction patterns are presented for selected material junctions and these 
are compared with data obtained by other means. Issues related to the numerical 
implementation are also discussed, including those pertaining to the convergence of 
the solution and sampling criteria. Finally, some family curves are given for selected 
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junction geometries in which the slab thickness is varied. 

Figures 5.5. 5.6, and 5.7 present the echowidth of several recessed slab geometries. 
These results were computed with the GSMF-GIBC solution of chapter three and the 
more general GIBC solution of chapters four and five. The slab having r = .04, e = 2. 

= 1.2 (see Figure 5.5) may be adequately simulated by a low contrast second order 
GIBC. whose solution is distinguished by the presence of a single unknown constant 
big. The computed value of this constant is given in Tables 5.2 and 5.3. which also 
contain the values of the constants pertaining the simulations given in Figures 5.5 
to 5.14. It is observed that the constants predicted by the GIBC and GIBC-GSMF 
solutions are practically identical, and the corresponding diffraction curves overlay 
one another for the PEC stub case as well as the PMC stub case. 

The solution constants for the GIBC solution corresponding to the PEC stub were 
determined by applying the boundary condition E y = 0 at the junction. In particular 
it was found that a satisfactory solution for bio could be obtained by enforcing this 
condition at a single point along the junction and for a single angle of incidence. 
The nearly-exact values in Table 5.2 result from enforcing the vanishing electric field 
at a single junction point for four distinct angles of incidence. We remark that 
there is no need for additional sampling points along the junction, since the interior 
$ function is constant with respect to y for a second order low contrast GIBC. A 
question may arrise, however, as to why it is desirable to sample at a greater number 
of sampling points than the number of unknowns. This is because the GIBC solution 
should ideally satisfy the boundary condition over all angles of incidence and at all 
points on the junction. Hence by using a sampling grid which spans the junction, 
it is possible to obtain a solution which satisfies the boundary conditions across 
the junction in an average sense. Once the overdetermined system is generated, a 


130 


solution for the constant(s) may be obtained by standard least-squares techniques. 

With regard to the computation of the field quantities at the edge, care must 
be exercised in evaluating the branch cut integrals. As seen from Table 5.1, the 
convergence of these integrals is not always guaranteed, since the integrand must 
have an asymptotic behavior of |A| -1-6 with 6 > 0. It may then be deduced from 
Table 5.1 that one cannot match E y or H y at the junction unless the order of all 
GIBC/GSTC for finitely-conducting bodies is even. We remark, however, that al- 
though the branch cut integrals for the PEC and PMC cases above always converge 
(they behave asymptotically as |A| and |A| respectively), their evaluation is 
not trivial by virtue of the infinite limits of integration. Herewith the infinite interval 
is transformed to a finite one. In addition, in case of a pole near the integration path, 
the addition and subtraction process described in [51] is employed to regularize the 
integrand. 

In Figure 5.6, the slab thickness is increased to .1A and it is now seen that 
fourth order conditions are required, resulting in three constants to be determined. 
It is lound that nearly 12 junction constraints are needed to adequately specify the 
constants by enforcing field continuity at three points across the junction for four 
angles of incidence. The agreement between the GIBC, GIBC-GSMF, and GSMF- 
exact are excellent, and the same is, of course, true for the GIBC and GIBC-GSMF 
constants (see Table 5.2). Note that the error in setting the constants to zero in this 
case is significant. The final recessed slab geometry has r = .4, t — 5, = 3 modeled 

by second order high contrast GIBCs. Again all solutions agree quite well as do the 
values of the constant 6 10 as given in Table 5.2. 

Figures 5.S to 5.14 present diffraction patterns for material half-planes of increas- 
ing thickness and are compared to data from a numerical model. The numerical 
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model was constructed by first generating the transient response of a finite length 
slab using bandlimited frequency domain data. The contribution from the half-plane 
edge was then extracted by time gating the transient response. This numerical model 
is valid except near grazing, where the surface wave and ray field excited by the back 
edge arrives at the front edge in concert with the incident plane wave. Figures 5.S to 
5.11 present the echowidth for low contrast simulations of a half-plane having e = 2, 
fi = 1.2. Clearly, the agreement between the numerical data and GSTC solution is 
excellent. The reader should also note the small values obtained for the constants as 
the thickness tends to zero, and this is in agreement with the second order GSTC so- 
lutions proposed in [53] and [54]. However, for thicker half-planes the constant plays 
a more significant role as evidenced by the erroneous result predicted in Figure 5.11 
when the constants are set to zero. 

Data based on two high contrast simulations are presented in Figures 5.12 and 
5.13 for a material half plane having e = 5 and n = 3. Because of the higher 
sampling required, numerical results could only be furnished for a thickness of up 
to .05A (see Figure 5.13). We observe that the results for the .01 A thick half-plane 
given in Figure 5.12 are in agreement and the same is generally true for the curves in 
Figure 5.13 despite the obvious instabilities of the numerical data. In Figure 5.14 a 
GSTC simulation is constructed for a two layer half- plane having T\ = .005, ei =5, 
/Xi = 3 and r 2 = .03, e 2 = 2, /x 2 = 1.2. The agreement of the GSTC solution with 
the numerical data is quite good except at edge on incidence, and judging by the 
abnormal behavior of the numerical data in this region it is conjectured that these 
data are in error. 

Finally, in Figures 5.16 to 5.21, family curves are given for various half-planes 
and grounded junctions. These are limited in thickness by computational restric- 
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tions stemming from the rapid increase in the number of constants as the order is 
increased, and additionally from the numerical intensiveness of the routine which 
determines the unknown constants. It is, therefore, necessary that the proper order 
(the one providing adequate simulation and converged results) of the GSTC be de- 
termined beforehand. This can be generally found by leaving out the constants and 
computing the diffraction coefficient for increasing order of GSTC until convergence 
w reached. The constants can then be determined for the order of the GSTC ren- 
dering convergence. This procedure was found quite adequate and was employed to 
generate the data in Figures 5.15 to 5.21. 
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Figure 5.5: //^polarization backscatter echo width for a recessed slab with r = . 4 , 

T ki ’/o7 1-2 modeled °( T ) second order low contrast GIBC (see 
table 5.2 for constants). v 
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Figure 5.6: H, polarization backscatter echo width for a recessed slab with r = .1, 
e = 2, fi = 1.2 modeled by 0(t 3 ) fourth order low contrast GIBC (see 
Table 5.2 for constants). 
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Figure 5.7: H z polarization backscatter echo width for a recessed slab with r = .4, 
t = 5, n = 3 modeled by 0(/c -1 ) second order high contrast GIBC (see 
Table 5.2 for constants). 
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Figure 5.8: ^polarization backscatter echo width for a material half-plane with 

T k. ' - o / = 2 ’ ^ = 1/2 modeled b y 0 { t ) low contrast GSTC (see 
lable 5.2 for constants). 
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Figure 5.9: H z polarization backscatter echo width for a material half-plane with 
t = .05, e = 2, fj. = 1.2 modeled by 0{r) low contrast GSTC (see 
Table 5.2 for constants). 
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Figure 5.10: H z polarization backscatter echo width for a material half-plane with 
r = .10, e = 2, fi = 1.2 modeled by O(r) low contrast GSTC (see 
Table 5.2 for constants). 
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Figure 5.11: H z polarization backscatter echo width for a material half-plane with 
r = .20, e = 2, n = 1.2 modeled by O(r) and 0(r 3 ) low contrast GSTC 
(see Table 5.2 for constants). 
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Figure 5.12: H z polarization backscatter echo width for a material half-plane with 
r = .01, t = 5, n = 3 modeled by 0 (ac°) (even symmetry) and 0(« -1 ) 
(odd symmetry) high contrast GSTC. The free space side is modeled 
by 0(t ) (even symmetry) 0(t 2 ) (odd symmetry) low contrast GSTC. 
(see Table 5.3 for constants). 
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Figure 5.13: H z polarization backscatter echo width for a material half-plane with 
t = .05, e = 5, n — 3 modeled by O(k 0 ) (even symmetry) and 0(k~ 1 ) 
(odd symmetry) high contrast GSTC. The free space side is modeled 
by O(t) (even symmetry) 0(r 2 ) (odd symmetry) low contrast GSTC. 
(see Table 5.3 for constants). 
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Figure 5.14: H z polarization backscatter echo width for a two layer .07A symmetric 
material half-plane with t x = .005, ei = 5, fi x = 3 and r 2 = .03, e = 2, 
M = 1*2. Layer 1 is modeled by 0(/c°) GSTC, layer 2 and the free space 
layer are modeled by O(t) (even symmetry) 0(r 2 ) (odd symmetry) low 
contrast GSTC (see Table 5.3 for constants). 
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Figure 5.15: H z polarization backscatter echo width curves for different thickness 
material half-planes with e = 2, fi = 1.2. All cases are modeled by O(r) 
low contrast GSTC (see Table 5.3 for constants). 
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Figure 5.16: H z polarization backscatter echo width curves for different thickness 
grounded slabs with e = 5, n = 3. The number to the left of the 
colon denotes the 0(t ) of the l.c. GIBC for the free space side, and the 
number to the right of the colon denotes the (9(«“ l ) of the h c slab 

GIBC. 
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Figure 5.17: H z polarization backscatter echo width curves for different thickness 
grounded slabs with e = 11, n = 7. The number to the left of the 
colon denotes the 0(t) of the l.c. GIBC for the free space side, and the 
number to the right of the colon denotes the 0(/c _1 ) of the h.c. GIBC 
for the right hand slab. 
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Figure 5.18: H : polarization backscatter echo width curves for different thickness 
grounded junctions, with ej = 2, fj.\ = 1.2 and £2 = 5, H 2 — 3. The 
number to the left of the colon denotes the O(r) of the l.c. GIBC for 
the left hand slab and the number to the right of the colon denotes the 
0(k~ 1 ) of the h.c. GIBC for the right hand slab. 
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Figure 5.19: H z polarization backscatter echo width curves for different thickness 
grounded junctions, with Ci = 5, = 3 and e 2 = 11, M 2 = 7. The 

number to the left of the colon denotes ue 0(«; _1 ) of the h.c. GIBC 
for the left hand slab and the number to the right of the colon denotes 
the 0(/c -1 ) of the h.c. GIBC for the right hand slab. 
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Figure 5.20: H z polarization backscatter echo width curves a grounded junctions, 
with ex = 2, fix = 1.2, Tj = .025 and e 2 = 5, = 3. The first number 

to the left of the colon denotes the O(t) of the l.c. GIBC for the left 
hand slab and the first number to the right of the colon denotes the 
0 (k *) of the h.c. GIBC for the right hand slab. The second number 
on either side of the colon denotes the O(r) of the l.c. GIBC for the 
free space slab needed to give the two sides the same thickness. 
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Figure 5.21: H x polarization backscatter echo width curves a grounded junctions. 

with ti = 5, Hi = 3, Tj = .05 and — 11, Hi = 7. The first number 
to the left of the colon denotes the 0(/c -1 ) of the h.c. GIBC for the 
left hand slab and the first number to the right of the colon denotes the 
0(k *) of the h.c. GIBC for the right hand slab. The second number 
on either side of the colon denotes the 0 ( t ) of the l.c. GIBC for the 
free space slab needed to give the two sides the same thickness. 
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FIGURE 5.5 
3IBC (GSMF) -PEC STUB 
tolO = -0 .0451 + j0.0212 
3I3C (GSMF) -PMC STUB 
blO - 0 . 6875+jO .3242 


GIBC-PEC STUB 

blO - -0 . 0451+ jO .0213 
GIBC-PMC STUB 

blO = 0. 6875+jO. 3242 


FIGURE 5.6 
3IBC (GSMF) 

blO = -0 . 457 8 + j 0 . 2593 
b20 = -0 . 0401-j0 . 0388 
3IBC 

blO = -0 . 4596+jO .2599 
b20 - -0 . 0400- jO . 0388 


bll = -0 . 0408- jO . 0111 
bll =- -0.041Q-j0.0110 


FIGURE 5.7 

3IBC (GSMF) GIBC 

blO = -0 . 0831+ jO .0489 blO = -0 . 0837+ jO . 048 9 

FIGURE 5.8 

3IBC 

blO .OOOO+j.OOOO 

FIGURE 5 . 9 

GIBC 

blO - . 0022- j .0009 
FIGURE 5.10 

GIBC 

blO - . 0008- j . 2323 

FIGURE 5.11 
GIBC 2nd order 

blO - 0 . 0152- jO .0166 
GIBC 4th order 

alO - -0 . 0734+ jO . 0138 all- 0 . 0145- jO . 0010 
a20 - -0 .0022+j0 .0103 

blO - 0 . 0350- jO .0368 bll - 0 . 0146+ jO . 0097 bl2 = 0.0035-j0021 

b20 - -0 . 0013+ jO . 0565 b21 - -0 . 0004+ jO . 0125 

b30 - -0 . 0022- jO .0015 


Table 5.2: Values of solution constants for curves presented in Figures 5 to 12 
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FIGURE 5 . ll 

alO = -0 . 0012- jO . OOOO 

blO =* 0. 0019- jO. 0005 

FIGURE 5.13 
alO - -0 . 0798- jO . 0054 

blO - 0 . 0171- jO . 0239 bll = 0 . 0090+ jO . 0002 

b2 0 = -0 . 0019 + j0 . 0011 

FIGURE 5 .14 

alO = 0 . 0003-j0 . 0000 

blO = 0 . 0180- jO . 0102 bll =■= -0 . 0060 + jO . 0002 

b20 - 0.0002-j0 .0001 

FIGURE5 . 15 

C-.005 

blO - -.0000+j .0000 
t= . 05 

blO - . 0022- j .0009 

t-.lO 

blO - . 0008- j .2323 

t= . 20 

blO - 0 . 0152- jO .0166 


Table 5.3: Values of solution constants for curves presented in Figures 12 to 15 


152 


5.5 Summary 

In this chapter a method was proposed and demonstrated for determining the 
unknown constants. Specifically, an eigenfunction expansion was presented as a 
representation for both the interior and exterior fields. The solution of chapter four 
was then recast into this form, allowing the unknown constants to be determined 
by imposing additional continuity conditions across the junction. Various scattering 
patterns were given validating the derived diffraction solution for several materia! 
half-planes and junctions. 



CHAPTER VI 


SUMMARY 


The goal of this dissertation was to develop a plane wave diffraction coefficient 
for thick multilayered symmetric slab junctions using simulations based on the gen- 
eralized impedance boundary conditions and generalized sheet transition conditions. 
This task was accomplished in four chapters dealing with the derivation of the 
GIBC/GSTC, the formulation and formal solution of the plane wave diffraction by 
a class of slab junctions, and the subsequent resolution of the unknown constants 
which arise in these solutions. 

To model multilayered slab junctions, it was first necessary to develop GIBC 
and GSTC for multilayered planar slabs, and this was the subject of chapter two. 
Recurrence relations were initially developed to relate fields in the adjacent layers 
of a multilayered structure, and these were subsequently used to derive infinite or- 
der boundary/transition conditions, conveniently expressed in a matrix product with 
each matrix corresponding to a layer. Low and high contrast approximations were 
then introduced to approximate the individual elements of each layer-matrix leading 
to a finite order GIBC/GSTC for the multilayered slab. Since each individual layer 
in the slab was characterized by a separate matrix, a low or a high contrast approx- 
imation could be employed for each individual matrix as dictated by the refractive 
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index of the layer. The accuracy of the derived conditions was evaluated by com- 
paring the GIBC/GSTC reflection coefficients to their exact counterparts and design 
curves were given for various single layer geometries based on the maximum phasor 
error of the approximation. 

In chapter three, the plane wave diffraction by a recessed slab in a ground plane 
was formulated and solved via the GSMF in conjunction with the dual integral 
equation method. The motivation for studying this geometry stemmed from the 
availability of exact data for uniform slabs. It thus served as a reference for evaluating 
the accuracy of GIBC in junction simulations. Furthermore, the employed GIBC 
simulation resulted in a unique solution, thus bypassing the non-uniqueness issue 
associated with GIBC/GSTC simulations of more arbitrary material junctions. 

The diffraction by a material discontinuity in a thick dielectric/ferrite slab was 
considered in chapter four. The slab was modeled by a distributed current sheet 
obeying generalized sheet transition conditions (GSTCs). This representation was 
then used to develop dual integral equations in terms of even and odd unknown 
spectral functions, which were proportional to the sheet currents. The solution for the 
spectra paralleled standard procedure but resulted in expressions involving unknown 
constants, revealing the non-uniqueness of the GSTC. It was demonstrated that the 
unknown constant(s) could be determined explicitly for the recessed slab discussed in 
chapter three by comparison with the results therein and it was also shown that the 
obtained solution reduces to simpler known solutions, including that for combinations 
of resistive and conductive sheets junctions. 

One way of determining the unknown solution constants discussed in chapter four 
is to employ field continuity across the junction. This however requires knowledge 
of the interior fields and in an effort to determine them, an eigenfunction decompo- 
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sition was developed for arbitrary symmetric multilayered slabs valid in the exterior 
and interior regions. The solution given in chapter four was then recast in this form 
permitting the determination of the interior slab field via analytic continuation. The 
unknown constants were then resolved by applying field continuity accross the junc- 
tion. Specifically, a point matching scheme was proposed in which the junction was 
sampled across its width for different angles of incidence. The resulting unique solu- 
tion was then verified for a number of recessed slab geometries by comparison with 
the GSMF solution in chapter three, and also for thin material half-planes whose 
scattering patterns were computed numerically. It was observed that the constants 
are very small for thin slab structures, but become significant with increasing slab 
thickness. Finally, family curves were generated for a number of half-plane and junc- 
tion geometries, and it was found that the convergence of the GIBC/GSTC solutions 
can be estimated by leaving out the solution constants. 

Clearly, the most challenging part of this work was the determination of the 
unknown solution constants. The mere fact that this was possible proved that the 
GSTC/GIBC are useful for practical simulations. Certainly, the method used in 
chapter five for determining the constants could be employed or paralleled in other 
applications. However, as noted therein, the evaluation of the constants for higher 
order GSTC simulations becomes numerically intensive and it would, therefore, be 
desirable to find alternative means for accomplishing this. For example, instead 
of point matching one could explore the orthogonality of the expansion modes or 
perhaps use a more efficient evaluation of the integrals. 

In this work we explored one application of the GSTC/GIBC simulation, that 
of diffraction by multilayered material slab junctions. As can be expected there 
are numerous other applications where the GSTC/GIBC can permit analytical solu- 
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tions. Examples include the extension of this work to skew incidence; possible char- 
acterization of junctions other than vertical where this solution may be employed 
in conjunction with a numerical one; and the diffraction by material junctions on 
curved surfaces at normal and skew incidences. Also, the characterization of multiple 
diffraction effects among material junctions is a straightforward process following the 
method already employed in [13]. 
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APPENDIX A 

MULTIPLICATIVE SPLIT FUNCTIONS 


In this appendix we consider the splitting of 

Q (A 2 ) = U A (A 2 ) + Vl — X 2 U B (a 2 ) (A.l) 

as a product of two functions, one of which is free of poles, zeros and branch cuts in 
the upper half of the A plane and the other having the same properties in the lower 
half of the A plane. That is, we seek to write Q (A 2 ) in the form 

e(A 2 ) =a + (A)S_(A) (A. 2) 


where the superscript + and - indicate an upper or lower function, respectively. 
Noting that 



U A (A 2 ) 

ii 

[l - A 2 ] 

l n 

(A. 3) 


Ub (A 2 ) 

K b 

= 

n= 0 

'l - A 2 

n 

(A. 4) 

with N a ■ 

= N b or N a = Nb + 1 , we 

may rewrite Q (A 2 ) as 



S(A’) = 

N S 





= E-Sn[^ 

n= 0 

- A 2 " 

i 

(A. 5) 

where N s 

= Afax(2N A ,2Ng + 1) and S n = A n / 2 if n is even and S n 

= B [n - 1 )/ 2 if n 

is odd. However, since we seek a multiplicative splitting of (A. 5), a more convenient 
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form to represent Q (A 2 ) is 


e M = so n (t + 


(A. 6) 


in which 7 n ‘ denote the zeros of the polynomial ( — A)*. We immediately 

identify that each of the product terms in (A. 6) can be factored as 


; now 


\/\ — A* def , 


7 


= M+ (A; 7) A/_ (A; 7) 


where 


A' + (A;l/ 7 ) = 


vl — A 


A. 1 


(A.S) 


A/+ (A; 7) 

is the split function characteristic to the impedance half plane having a constant 
surface impedance I/7 [34]. With the branch choosen so that Im{\/l — A 2 ) < 0, 
M + (A; 7) is explicitly given by 

f M+ (A; 7) /m(7) < 0 


A/+ (A; 7) = M- ( — A; 7) = 


JV 




Im( 7) > 0, 


M + (cos a; I/77) 
In this. 


»n(*/2) [i + V2cos( ,, ey , ' l, )l [1 + Vicos 


75 [«, (3jt/ 2 -£>-«)*„ (x/2 - a + «)] ! 


(A.9) 


(A.10) 


Im (77) > 0 
A = cos a 

Im (y/l ~ l/ 7 ? 2 ^ — 0 
0 = sin _1 (77) with 0 < Re(9 ), 

(A. 11) 


and v f I (a) is the Maliuzhinets function [22] whose evaluation in algebraic form has 
been given in [50]. We remark that in the limit as 7 — ► 0, 

yj T^A 2 


M+ (A; 7 *• 0) = 


Vi 


(A. 12) 
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and as 7 — ► 00 

M+ (A; 7 — ► oc) = 1 (A. 13) 

The determination of ( A) is now rather trivial. By substituting (A. 7) into 
(A. 6) we easily obtain 

e + (A) = r.(-A) = (a.u> 

n= 1 
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